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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (SI)  units  as  follows: 


Multipl 


cubic  feet  per  second 


miles  (U.  S.  statute) 


0.02831685 

0.3048 

1.609344 


To  Obtain 


cubic  metres  per  second 

metres 

kilometres 
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DEVELOPMENT  OF  A  NUMERICAL  MODELING  CAPABILITY  FOR  THE 
COMPUTATION  OF  UNSTEADY  FLOW  ON  THE  OHIO  RIVER 
AND  ITS  MAJOR  TRIBUTARIES 


PART  I :  INTRODUCTION 
Purpose 

1.  The  U.  S.  Army  Engineer  Division,  Ohio  River  (ORD),  is  respon¬ 
sible  for  the  multipurpose  operation  of  many  reservoirs  in  the  Ohio 
River  Basin.  During  periods  of  flooding  on  the  lower  Ohio  and  lower 
Mississippi  Rivers,  the  Division  Office  issues  specific  instructions  for 
the  operation  of  Barkley  and  Kentucky  Reservoirs  on  the  Cumberland  and 
Tennessee  Rivers,  respectively.  The  primary  objective  of  flood-control 
regulation  by  such  reservoirs  is  to  reduce  the  frequency  and  magnitude 
of  flooding  of  lands  along  the  Ohio  River  and  its  tributaries.  In 
addifion,  ORD  also  maintains  and  operates  navigation  locks  and  dams  on 
the  Ohio  River.  At  the  present  time,  these  structures  are  operated  only 
for  navigation  purposes;  however,  studies  are  being  considered  to  assess 
the  feasibility  of  utilizing  the  pools  for  hydropower.  Obviously,  a 
mathematical  model  capable  of  computing  unsteady  flows  on  the  Ohio  River 
and  its  major  tributaries  for  different  tributary  inputs  while  allowing 
for  the  influence  of  the  operation  of  the  navigation  locks  and  dams  on 
the  Ohio  River  would  be  a  useful  tool  to  aid  in  these  activities.  The 
development  of  such  a  modeling  capability  for  the  ORD  was  the  objective 
of  this  study. 


Background 

2.  The  study  reported  herein  was  in  essence  an  extension  of  a 
modeling  study  that  began  in  1971.  At  its  twenty- seventh  meeting  on 
19  May  1970,  the  Mississippi  Basin  Model  (MBM)  Board  authorized  a  study 
to  develop  computer  programs  for  unsteady  flow  computations  along 
reaches  of  the  Mississippi  River  and  its  larger  tributaries.  At  the 
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thirty-second  meeting  on  7  January  1971,  in  a  joint  effort  with  ORD, 
the  U.  S.  Army  Engineer  Waterways  Experiment  Station's  area  of  re¬ 
sponsibility  was  determined  to  be  the  lower  Ohio  River  from  Louisville, 
Kentucky,  to  the  junction  with  the  Mississippi  River  at  Cairo,  Illinois. 
Results  from  this  effort  were  reported  by  Johnson  and  Senter  (1973)  and 
illustrated  the  applicability  of  a  finite-difference  representation  of 
the  unsteady  flow  equations  applied  to  the  lower  Ohio  River. 

3.  At  the  conclusion  of  the  study  by  Johnson  and  Senter,  it 
was  decided  that  in  order  to  provide  ORD  with  a  model  capable  of  being 
used  to  predict  stages  at  Cairo  for  given  release  schedules  at  Barkley 
and  Kentucky  Reservoirs,  the  Cumberland  and  Tennessee  Rivers  should  be 
treated  as  dynamic  branches  of  the  system  (see  Figure  1  for  a  schematic 
of  the  modeled  region).  In  the  earlier  study,  these  flows  were  treated 
as  lateral  inflows  into  the  Ohio;  therefore,  a  study  was  initiated  to 
improve  the  modeling  capability  to  allow  for  the  modeling  of  flows  in 
multijunction  systems.  In  addition,  it  was  also  considered  necessary  to 
obtain  more  accurate  geometric  data  from  the  MBM.  The  one-dimensional 
mathematical  mo^el  that  resulted  from  this  study  was  an  extension  of  a 
model  developed  at  the  Tennessee  Valley  Authority  (TVA)  by  Garrison, 
Granju,  and  Price  (1969).  A  brief  discussion  of  the  model  (called 
SOCHMJ  from  Simulation  of  Open  Channel  Hydraulics  in  Multi-Junction 
Systems)  and  results  from  the  study  are  presented  by  Johnson  (1974). 

4.  As  a  result  of  the  success  enjoyed  from  the  use  of  SOCHMJ  as 
an  aid  in  determining  the  operation  of  Barkley  and  Kentucky  Reservoirs 
during  periods  of  flooding  on  the  lower  Ohio  and  Mississippi  Rivers,  ORD 
then  initiated  a  project  to  extend  the  modeling  limits  to  McAlpine  Lock 
and  Dam  on  the  Ohio,  with  the  Green  and  Wabash  Rivers  being  treated  as 
dynamic  branches  of  the  system.  This  effort  required  developing  a 
capability  in  SOCHMJ  for  handling  the  high-lift  locks  and  dams  at  Can- 
nelton,  Newburgh,  and  Uniontown  as  illustrated  in  Figure  2.  Results  and 
a  discussion  of  the  lock  and  dam  routine  as  well  as  other  modifications 
to  SOCHMJ  are  presented  by  Johnson  (1977). 

5.  The  extension  of  SOCHMJ  to  McAlpine  Lock  and  Dam  resulted  in  a 


J 
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model  which  simulated  approximately  750  river  miles, *  and  computation 
costs  suggested  that  this  size  system  approached  the  economic  limit  for 
the  numerical  scheme  used  in  SOCHMJ.  This  stems  from  the  fact  that 
SOCHMJ  utilizes  an  explicit  finite-difference  scheme  to  solve  the 
governing  equations.  As  a  result,  the  computational  time-step  At  is 
restricted  by  the  stability  criterion 


where  all  variables  are  defined  in  Appendix  C.  For  values  of  Ax  of  3 
to  5  miles  on  the  Ohio  River,  the  time-step  is  restricted  to  a  maximum 
of  perhaps  5  min.  With  such  a  small  computational  step,  computational 
costs  for  a  larger  system  would  become  excessive. 

6.  After  the  extension  of  SOCHMJ  to  McAlpine  Lock  and  Dam,  the 
ORD  expressed  an  interest  in  developing  a  modeling  capability  for  the 
complete  Ohio  River  from  Pittsburgh,  Pennsylvania,  through  its  junction 
with  the  Mississippi  River.  Based  upon  results  from  the  previous  study, 
it  was  obvious  that  the  use  of  SOCHMJ  was  not  economically  feasible. 
Therefore,  the  development  of  a  new  modeling  capability  that  would  not 
contain  such  a  restrictive  stability  criterion  was  required. 

Scope 

7.  When  employing  the  method  of  finite  differences  to  solve  the 
governing  equations  there  are  two  general  types  of  schemes,  either  an  ex¬ 
plicit  or  an  implicit  scheme.  As  previously  noted,  explicit  schemes  are 
conditionally  stable,  whereas  implicit  schemes  are  not.  In  other  words, 
much  larger  computational  time-steps  may  be  employed  with  an  implicit 
scheme,  which  of  course  in  general  reduces  the  cost  of  computations. 

8.  A  one-dimensional,  implicit  finite-difference  model  developed 
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*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  4. 
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by  Chen  (1973)  and  applied  to  a  segment  of  the  upper  Mississippi  River 
by  Simons  et  al.  (1975)  was  selected  as  the  base  from  which  an  efficient 
modeling  capability  for  the  complete  Ohio  River  Basin  (Figure  3)  could 
be  developed  and  provided  to  ORD.  Necessary  modifications  included  the 
development  of  input/output  routines,  the  generalization  of  the  basic 
code  to  handle  a  system  containing  an  unlimited  number  of  branches  and 
junctions,  and  the  incorporation  of  a  technique  similar  to  that  pro¬ 
grammed  in  SOCHMJ  to  handle  the  many  high-lift  locks  and  dams  in  the 
system.  The  resulting  computer  code  is  called  FLOWSED  to  reflect  the 
fact  that  both  sediment  movement  and  flow  are  computed.  The  work  re¬ 
ported  herein  emphasizes  the  flow  computations  with  coefficients  selected 
to  make  the  sediment  transport  negligible. 

9.  In  addition  to  the  work  required  on  the  basic  computer  program 
furnished  by  Chen  (1973)  to  convert  it  into  an  efficient  model  for  use 
on  the  Ohio  River  and  its  major  tributaries,  necessary  geometry  data  had 
to  be  constructed.  These  data  were  constructed  from  storage  volume  data 
collected  from  the  MBM,  topographic  maps  and  cross-sectional  information 
furnished  by  ORD.  To  convert  the  data  obtained  from  the  latter  two 
sources  into  the  form  required  by  FLOWSED,  a  computer  program  called 
GEOM  was  written. 

10.  After  the  generation  of  the  required  tables  of  geometry  data, 
historical  flood  data  from  1964,  19/2,  and  1976  were  used  to  calibrate 
various  portions  of  the  system.  Such  a  calibration  requires  the  manipu¬ 
lation  of  Manning's  n  and/or  geometry  data  until  recorded  and  computed 
elevations  and  discharges  are  in  reasonable  agreement. 

11.  PART  II  of  this  report  presents  a  discussion  of  the  governing 
equations  to  be  solved  and  underlying  assumptions  made  in  their  deriva¬ 
tion,  along  with  the  particular  solution  technique  developed  by  Chen 
(1973).  In  addition,  special  features  such  as  the  handling  of  naviga¬ 
tion  dams  and  levee  overtopping  are  also  discussed.  A  general  descrip¬ 
tion  of  the  computer  code  with  special  emphasis  on  the  computational 
cycle  and  input  requirements  is  presented  in  PART  III.  In  PART  IV, 
details  of  the  complete  Ohio  River  system  and  construction  of  the  re¬ 
quired  geometry  data  base  are  presented.  PART  V  presents  results  from 


the  calibration  effort,  while  PART  VI  summarizes  the  modeling  effort  and 
offers  recommendations  for  additional  work  needed.  A  listing  of  the 
input  cards  required  by  FLOWSED  is  given  in  Appendix  A;  Appendix  B 
presents  an  example  of  FLOWSED  output. 


PART  II:  THEORETICAL  DISCUSSIONS 


Governing  Equations 


12.  In  the  case  of  one-dimensional  open-channel  flows  within 
rigid  boundaries,  the  flow  behavior  can  be  adequately  described  by 
the  Saint  Venant  partial  differential  equations  of  unsteady  flow. 

These  equations  are  derived  by  considering  the  conservation  principle 
for  mass  and  for  the  momentum  of  the  flows.  For  mobile  boundary  chan¬ 
nels,  three  equations  are  needed  to  describe  the  sediment-water 
mixture:  continuity  and  momentum  of  the  mixture  and  conservation  of 
mass  of  the  sediment.  A  detailed  derivation  of  these  equations  is 
presented  by  Chen  (1973).  Some  of  the  more  important  features  are 
discussed  below. 

13.  The  governing  equations  are  derived  for  nonprismatic  channels 
with  irregular  cross  sections.  The  flow  field  under  consideration  is 
assumed  to  be  one -dimensional  and  the  pressure  field  varies  in  the 
vertical  direction  in  a  hydrostatic  manner.  These  assumptions  imply 
that  the  river  reach  to  be  modeled  should  be  reasonably  straight  with 
the  free  surface  taken  to  be  a  horizontal  line  across  the  section  and 
that  vertical  accelerations  are  negligible.  In  addition,  the  density  of 
the  water- sediment  mixture  is  homogeneous  in  the  vertical  direction. 
Sediment  continuity  equation 

14.  Considering  Figures  4a  and  4b,  let  Q  be  the  total  volume 

s 

of  sediment  transported  by  the  riverflow  per  unit  time,  where  Qg  is  a 
function  of  both  the  longitudinal  distance,  x  ,  and  time,  t  Let  q 

s 

be  the  unit  rate  of  volume  of  sediment  entering  the  river  because  of 

lateral  inflows  into  the  river  as  a  result  of  small  tributaries  or 

overland  flow,  where  q  is  expressed  in  volume  of  sediment  per  unit 

s 

length  and  unit  time.  Equating  the  change  in  the  mass  of  sediment 
stored  within  the  control  volume  presented  in  Figure  4b  during  the  time 
interval  At  to  the  difference  between  the  sediment  entering  and  the 
sediment  leaving  results  in  the  sediment  continuity  equation,  a  state¬ 
ment  of  which  is 
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Sediment  entering  -  Sediment  leaving  =  Change  in  sediment  stored  (1) 


15.  The  change  in  the  storage  of  sediment  within  the  control 
volume  is  effected  in  two  ways;  by  the  deposition  or  scour  of  sediment 
on  the  riverbed  and  by  the  change  in  the  suspended  sediment  concentra¬ 
tion.  Assuming  that  the  deposition  or  scour  occurs  uniformly  over  the 
whole  bed  area,  the  mass  of  sediment  per  computational  cell  in  a  bed 

layer  of  thickness  A2  is  p  (PAZAx)p  where  p  is  density  of 

s  s 

sediment,  P  is  the  wetted  perimeter  at  the  section  where  the  control 
volume  is  located,  and  p  is  the  porosity  of  the  bed  layer,  i.e.  the 
volume  of  sediment  per  unit  volume  of  the  bed  layer.  The  change  in 
sediment  storage  due  to  the  change  in  the  suspended  sediment  concentra¬ 
tion  taken  as  an  average  over  the  control  volume  can  be  expressed  as 

p  (9/9t)(AAxC  )At  ,  where  A  is  the  cross-sectional  flow  area  and  C 
s  s  s 

is  the  suspended  sediment  concentration. 

16.  Assuming  that  AZ  =  (9Z/3t)At  ,  the  sediment  continuity  equa¬ 
tion  expressed  as  Equation  1  becomes 


9Q  ~7  3AC 
s  vZ  ,  s 

9x~  +  P  P  + 


9t 


= 


where  Z  is  the  elevation  of  the  channel  bed. 


(2) 


Continuity  equation  for 
sediment-water  mixture 

17.  In  this  case,  both  the  mass  of  water  and  sediment  are  con¬ 
sidered  together.  As  for  the  sediment  continuity,  the  basic  equation 
is  expressed  over  the  control  volume  as 


(Sediment-water  mixture  entering)  -  (Sediment-water  mixture  leaving) 

=  Change  in  storage  of  sediment-water  mixture  (3) 


It  can  be  shown  that  if  Q  is  the  total  discharge,  A  is  the  total 
cross-sectional  area,  and  q^  is  the  unit  lateral  inflow  and  if  Equa¬ 
tion  2  is  utilized,  then  the  mathematical  representation  of  Equation  3 
becomes 


(4) 


r  ■ 


t 

3 


as 

dx 


+ 


§a  p  az 
at  v  at 


=  «i 


Momentum  equation  for 
sediment-water  mixture 

18.  For  the  control  volume  in  Figure  4a,  the  equation  for  the 
conservation  of  momentum  of  the  sediment-water  mixture  can  be  written  as 

Net  rate  of  momentum  flux  into  the  control  volume  +  Sum  of  forces 

=  Rate  of  accumulation  of  momentum  within  the  control  volume  (5) 

Once  again  referring  to  Chen,  the  above  equation  can  be  written  in 
mathematical  terms  as 

§£  (pQ)  ^x  =  -  (pp  Ax  -  pq£  |  Ax  +  J  forces  (6) 

where  the  forces  acting  on  the  control  volume  are  gravity,  pressure,  and 

frictional  resistance,  p  is  the  density  of  the  sediment-water  mixture, 
and  p  is  the  momentum  correction  factor. 

Gravity  force 

19.  The  force  due  to  gravity  is  the  weight  of  the  fluid  mixture 

within  the  control  volume.  If  Sx  is  the  slope  of  the  control  volume 

bottom,  then  the  component  of  this  force  along  the  flow  direction 
becomes 


pgAAxSx  (7) 

if  the  flow  within  the  control  volume  is  uniform. 

Pressure  force 

20.  The  net  pressure  force  along  the  direction  of  flow  is  com¬ 
posed  of  the  pressure  forces  acting  on  the  two  ends  of  the  control  volume 
and  the  pressure  forces  in  the  direction  of  the  flow  on  the  banks  due 
to  widening  or  narrowing  along  the  length  of  the  nonprismatic  channels. 
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Assuaing  a  hydrostatic  pressure  distribution,  it  can  be  shown  that  the  net 
pressure  force  acting  on  the  control  volume  is 


-pgA  gj  Ax  (8) 

where  Y  is  the  flow  depth. 

Frictional  force 

21.  The  frictional  force  resisting  the  motion  of  the  sediment- 
water  mixture  acts  along  the  solid  boundaries  of  the  channel  and  is 
expressed  as  -(tqPAx)  ,  where  is  the  shear  stress  and  P  is  the 

wetted  perimeter.  Assuming  that  the  expression  for  the  shear  stress  in 
steady  flow  applies  in  unsteady  flows,  the  frictional  force  along  the 
direction  of  flow  acting  on  the  control  volume  can  be  expressed  as 

-pgASjAx  (9) 

where  is  the  friction  slope.  An  eddy  loss  term  Sg  similar  to  the 

friction  slope  has  been  added  to  Chen's  basic  equation  to  account 

for  head  losses  in  addition  to  those  due  to  boundary  friction.  These 
losses  are  due  to  large-scale  eddies  formed  in  the  flow  at  rather  abrupt 
changes  in  the  cross  sections  along  the  channel.  The  eddy  loss  slope  is 
evaluated  using  an  expression  taken  from  Fread  (1976): 


S 

e 


2g 


(10) 


in  which  K£  is  a  nondimens ional  coefficient  of  contraction  and  expan¬ 
sion.  Cross  sections  contracting  abruptly  in  the  direction  of  flow  have 
Kg  values  ranging  from  0  to  0.4,  while  abruptly  expanding  cross  sections 
have  Kg  values  ranging  from  0.5  to  1.0.  The  larger  coefficients  are 
associated  with  the  more  abrupt  contractions  and  expansions  of  the  cross 
section. 

22.  Considering  Figure  4a,  the  flow  cross-sectional  area  is  given 


Using  the  Leibnitz  rule,  one  can  then  compute  the  time  and  spatial  deriva 
tive  of  A  in  Equations  2,  4,  and  6  as 

■.♦tC.ng-.g  <«> 

o 

and 

Y 

3x  =  /  5i  5  *  J(x>r)  5J  =  4x  *  B  51  (13) 

O 

Y 

where  B  is  the  top  width  of  the  channel  and  A^  is  the  rate  of  change 
of  area  with  respect  to  x  with  the  flow  depth  Y  held  constant. 

23.  Substituting  the  expressions  for  the  gravity,  pressure,  and 
frictional  forces  as  well  as  the  expression  for  3A/3x  into  Equation  6, 
the  conservation  of  momentum  equation  for  the  sediment-water  mixture 
becomes 

It  <«>  ♦ v  s  +  f*Pv  i  -  W2*  s 

♦  *A  |j  (pY)  =  pgA(Sx  -  Sf  -  S.)  ♦  (SpvV  (14) 

24.  Equations  2,  4,  and  14  are  the  equations  governing  the  motion 
of  a  sediment -water  mixture  in  open  channels  in  a  one-dimensional  sense. 
These  equations  involve  six  unknowns;  namely,  the  flow  discharge  Q  ,  the 
flow  depth  Y  ,  the  bed  elevation  Z  ,  the  sediment  transport  rate  Qg  , 
the  frictional  slope  ,  and  the  density  of  the  sediment-water  mix¬ 
ture,  p  .  Other  variables  such  as  the  lateral  inflows,  bed  porosity, 


and  geometry  data  are  expected  to  be  known.  To  achieve  closure  of  the 
system,  three  additional  relations  are  required.  These  are  provided  by 
a  sediment  transport  function  given  as 


Qs  =  CjV  Y  JQ  (15) 

Manning's  equation  which  relates  the  friction  slope  to  the  flow  and  chan¬ 
nel  characteristics,  and  an  equation  of  state  for  the  density  of  the 
sediment-water  mixture  given  as 

p=p  +  C  (p  -p)  (16) 

v  Kw  s  rs  rw 

With  these  additional  relations  one  can  then  solve  for  the  basic  unknowns 
Q  ,  Y  ,  and  Z  . 


Special  Features 


Locks  and  dams 

25.  To  account  for  the  effect  of  navigation  locks  and  dams,  the 
following  equations  are  utilized  to  simulate  the  sediment-water  mixture 
through  the  locks  and  dams: 


Qs 

us 

=  Qs . 
ds 

(17) 

Qus 

=  Qds 

(18) 

* 

=  f(t) 

(19) 

where  the  subscripts  us  and  ds  refer  to  the  cells  immediately  up¬ 
stream  and  downstream  of  the  structure.  The  normal  procedure  is  to 
input  a  constant  water-level  elevation  upstream  of  a  lock  and  dam  to 
reflect  the  pool  elevation  the  lock  operator  is  expected  to  maintain, 
i.e.  f(t)  =  constant  .  With  such  a  procedure,  the  flow  required  to  be 
passed  through  the  structure  in  order  to  maintain  the  upstream  elevation 


is  computed.  Theoretically,  the  operator  could  then  use  gate  rating 
curves  to  make  the  gate  adjustments  required  to  pass  the  computed  flow. 
In  hydropower  feasibility  studies,  one  may  wish  to  specify  some  time 
variation  of  a  particular  pool  rather  than  prescribing  a  constant  value. 
Junctions 


26.  The  interaction  between  the  main  river  and  a  tributary  being 
handled  in  a  dynamic  fashion  as  opposed  to  being  treated  as  lateral  in¬ 
flow  is  simulated  by  the  following  continuity  and  energy  equations: 


q3  =q2  +  Q! 


(20] 


Qs  =QS  +  % 

S3  s2  Sl 


(21] 


+  Y2  +  2i  =  Z3  +  Y3  +  “2  4  +  (SfAx)2 


2g 


(22] 


2  2 
V  V 

Z,  +  Y,  +  ^  =  Z,  +  Y„  +  0,  ^  +  (SfAx)] 


'1  *1  2g  3  3  “1  2g 


(23] 


in  which  a  is  the  energy  correction  factor,  SfAx  is  the  energy  head 
loss  and  subscripts  1,  2,  3,  as  illustrated  below,  refer  to  sections 
above  the  confluence  on  the  tributary  and  main  river  and  below  the  con- 


Levees 


27.  To  handle  the  effect  of  levee  overtopping,  the  weir  equation 


q  = 


(24) 


m 


where  H  is  the  height  of  water  over  the  levee,  is  invoked  to  compute 
the  unit  lateral  outflow  from  the  channel  whenever  the  water-surface 
elevation  exceeds  the  height  of  the  levee.  A  basic  assumption  in  the 
current  handling  of  levees  is  that  flow  which  leaves  the  channel  as  a 
result  of  levee  overtopping  is  lost  from  the  system. 

Overbanks 

28.  The  lateral  inflow  q^  consists  of  two  components,  q^  and 

q.  ,  induced  by  the  handling  of  floodplains  and  tributaries,  respec- 
x’2 

tively.  For  overbank  flow  on  floodplains,  q.  is  computed  from 


AxAt 


Ah 


(25) 


where  Af  is  the  surface  area  of  the  floodplain  and  Ah  is  the  change 
over  a  time  period  At  of  the  water-surface  elevation.  Similarly,  the 
unit  lateral  inflow  of  sediment  is  also  composed  of  two  parts.  As  the 
water  surface  rises,  flow  onto  the  floodplain  carries  sediment  along 
with  it.  The  coarse  material  is  deposited  along  the  riverbank,  re¬ 
sulting  in  an  increase  of  the  height  of  the  natural  levee.  Chen  assumes 
a  triangular-shaped  natural  levee  with  bottom  angles  of  30  (face  to  the 
main  channel)  and  15  deg.  The  increase  in  height  of  the  natural  levee 
over  a  time  At  can  then  be  computed  from 


(26) 


where  the  sediment  concentration  near  the  riverbank  is  assumed  to  be 

five  times  larger  than  the  average  concentration  C  .  Therefore,  the 

s 


•  •  • 
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lateral  sand  flow  associated  with  overbank  flow  is  determined  from 


%  =  S’*,0, 


(27) 


on  the  rising  side  of  the  hydrograph,  but  is  set  to  zero  on  the  falling 
limb  as  water  returns  to  the  channel. 

29.  In  addition  to  the  movement  of  sandy  material,  the  stream  also 
carries  a  wash  load  composed  of  mainly  clays  and  silts.  Chen  has  selected 
the  following  expression  for  the  wash-load  discharge  on  the  upper 
Mississippi  River, 


=  3.53  x  10  iJQ 


13rt2. 752 


(28) 


The  deposition  of  silts  and  clays  on  the  floodplain  is  then  computed  from 

2(zytc wl)^ 


(29) 


where  -  QyjVQ  for  q^  <  0  but  is  set  to  zero  when  q^  >  0  .  Ad- 

ditional  discussion  of  the  manner  in  which  sediment  deposition  on  over¬ 
bank  areas  is  handled  is  provided  by  Simons  et  al.  (1975). 

Boundary  Conditions 

30.  The  governing  Equations  2,  4,  and  14  constitute  a  set  of 
first  order  hyperbolic  partial  differential  equations.  As  long  as  the 
flow  regime  is  subcritical,  i.e. 


V  <  yJ~gY 


(30) 


the  following  boundary  conditions  are  possible  for  a  well-posed  problem. 
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At  an  upstream  open  boundary,  either  the  flow  discharge  or  the  water- 
surface  elevation  may  be  prescribed  as  a  function  of  time.  In  addition, 
the  sediment  discharge  hydrograph  must  either  be  specified  or  can  be 
computed  internally  from  the  given  transport  function.  At  a  downstream 
boundary,  either  the  flow  discharge,  water-surface  elevation,  or  a  rating 
curve  relating  the  flow  discharge  and  elevation  must  be  prescribed. 

Unless  the  downstream  boundary  coincides  with  a  control  structure,  the 
normal  procedure  is  to  employ  a  rating  curve  as  the  downstream  boundary 
condition. 

Initial  Conditions 

31.  As  previously  noted,  the  equations  to  be  solved  constitute 
a  hyperbolic  system  and  initial  values  of  the  dependent  variables  must 
be  prescribed  to  begin  the  time  marching  of  the  solution.  A  steady- 
flow  profile  or  perhaps  a  transient  profile  from  previous  computations 
can  be  used.  The  specification  of  initial  conditions  is  flexible  due 
to  the  characteristic  of  hyperbolic  equations  that  the  solution  be¬ 
comes  independent  of  initial  conditions  after  a  sufficient  length  of 
time. 

Finite-Difference  Representation 


32.  The  governing  equations  expressing  the  conservation  of  sedi¬ 
ment  and  the  conservation  of  mass  and  momentum  of  the  sediment-water 
mixture  do  not  in  general  possess  analytic  solutions.  One  must  there¬ 
fore  rely  upon  numerical  techniques  such  as  finite  differences  to  obtain 
solutions.  The  finite-difference  approximations  selected  by  Chen  to 
express  the  values  and  the  partial  derivatives  of  a  function  4>  ,  where 
0  represents  the  dependent  variables  Q  ,  Y  ,  Z  ,  are  given  as 


4»n  , 

vi  +  vi+l 


9$  a 

dx  “ 


.n+1  ,n+l 

♦iti  :  h 

Ax 


(32) 


S 


1 


3$ 

at 


■it  [(♦r  ■  <)  ■  ♦•«)] 


(33) 


Figure  5  illustrates  the  positioning  of  the  derivatives  in  a  computa¬ 
tional  cell. 

33.  Constructing  difference  equations  from  the  governing  differ¬ 
ential  equations  through  use  of  the  finite-difference  approximations 
above  results  in  a  linear-implicit  finite-difference  scheme  in  which 
the  difference  form  of  the  governing  equations  written  over  a  computa¬ 
tional  cell  formed  by  sections  i  and  i+1  take  the  form 

„n  qQ+1  +  „n  „n+l  +  vn  Qn+1  +  if11  vn+^  +  v11  7n+l  _  pQ  /-oa\ 

KklQi  +  Kk2  i  +  Kk3Zi  +  Kk4Qi+l  +  Vi+1  +  Vi+1  Ek  t34) 

K?iQf 1  +  k“,y"+1  +  K^zf1  +  K“AQ“1!  +  K^cY“t;  +  =  E“  (35) 


22i 


23  i 


£4xi+l  25  i+1  £6  i+1 


KnlQn+1  +  Kn„Yn+1  +  K%Zn+1  +  K\Qn:;  +  K“  Y"!}  +  Kn,Z?!l  =  En  (36) 
mlxi  m2  i  m3  x  n»4xi+l  m5  l+l  m6  l+l  m 


where  k  =  3i  ,  £  =  3i+l  ,  m  =  3i+2  ,  and  the  coefficients  K's  and 

E's  are  functions  of  known  variables  from  the  previous  time  line. 

34.  In  casting  the  difference  equations  into  the  form  shown,  it 

might  be  noted  that  the  sediment  transport  function  Q  and  the  fric- 

s 

tion  slope  are  expanded  in  a  first  order  Taylor  series,  i.e. 

C  *  <C  ♦  ©V1  -  «”)  ♦  (§r)°  ■  *■)  ♦(£)"  (*an  -  z“) 


Sf1  3  S 


\  ♦  (  £)V‘  -  0")  +(§r)n  -  v")  ♦  (£*)“  U"tl  -  z”) 
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4  4 


in 


o 


where  the  partial  derivatives  are  computed  by  differentiating  the 
analytic  expressions  for  Qg  and  .  The  friction  slope  is  taken  at 
the  (n+1)  time  level  in  the  momentum  equation  to  assure  the  stability 
of  the  scheme. 

35.  Before  briefly  discussing  the  solution  technique  employed  by 
Chen  for  solving  the  algebraic  system  of  equations  resulting  from  Equa¬ 
tions  34-36,  the  manner  in  which  boundary  conditions  are  handled  is 
discussed.  Assume  that  the  flow  discharge  is  specified  at  an  upstream 
boundary,  i.e. 


Q.  =  f(t) 

this  equation  can  be  cast  into  the  form 


u 


K“Qn+1  +  K"  Yn+1  +  K?,Zn+1  =  E?+1 
k4xi  k5  1  k6  i  k 


(37) 


where  =  1  ,  K^,.  =  =  0  ,  and  E^  =  f(t)  .  The  subscript  k 

has  the  value  k  =  3 ( i- 1 )  +  1  and  thus  will  be  equal  to  1  at  the 
first  upstream  section.  As  can  be  seen,  Equation  37  is  in  a  form  similar 
to  the  difference  Equations  34-36.  Thus,  it  should  be  obvious  that 
boundary  conditions  can  be  set  as  part  of  the  overall  solution  scheme 
merely  by  setting  the  proper  coefficients  equal  to  either  one,  zero,  or 
the  proper  time  varying  function  on  the  right-hand  side. 

36.  When  specifying  a  sediment  discharge  hydrograph  at  an  up¬ 
stream  boundary,  i.e. 


V  * 


once  again  a  Taylor  series  expansion  is  utilized  to  cast  the  equation 
into  the  form  of  the  difference  Equations  34-36.  Therefore, 


^  \9Y  )i  ' 


yll+l 

i 
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is  rearranged  as 


Vi  +  Vi  +  Vi  =  E£ 


in  which  the  subscript  i  takes  on  its  value  at  the  upstream  boundary  and 
£  =  3(i-l)  +  2  . 

37.  At  a  downstream  boundary  where  a  rating  curve  is  prescribed, 
a  Taylor  series  expansion  once  again  enables  one  to  obtain  the  proper  form 
of  the  equation.  Assuming 


«i  =  W 


where  i  takes  its  value  at  the  downstream  boundary  and  lu  = 
is  the  water-surface  elevation,  one  can  write 


*  *  o?  *  (§r).(Yf 1  • Y?)  t(sr).(zf 1  -  z°i) 


which  can  be  rearranged  as 


K  ,Q.  +  K  _Y.  +  K  -Z.  =  E 
mil  m2  i  m3  i  m 


with  m  =  3i  .  If  the  elevation  rather  than  a  rating  curve  is  specified, 

one  would  set  K  ,  =  0  and  K  „  =  K  „  =  1 '  .with  E  being  set  to  the 
ml  m2  m3  m 

particular  time  variation  specified. 


Solution  Technique 

38.  The  solution  technique  developed  by  Chen  is  what  is  commonly 
called  the  double  sweep  algorithm.  To  demonstrate  the  procedure,  consider 
the  following  problem 


•  •••••••• 


composed  of  a  main  river  and  one  tributary.  The  set  of  linear  equations 
for  branch  1  becomes 


+  K1,5Y1  +  K1,6Z1  =  E1 


*2,4*!  +  K2,5Y1  +  K2,6Z1  =  E2 


K3,1Q1  +  K3,2Y1  +  K3,3Z1  +  K3,4Q2  +  K3,5Y2  +  K3,6Z2  =  E3 


K4,1Q1  +  K4,2Y1  +  K4,3Z1  +  *4,4*2  +  K4,5Y2  +  K4,6Z2  =  E4 


K5,1Q1  +  K5,2Y1  +  K5,3Z1  +  *5,4*2  +  K5,5Y2  +  K5,6Z2  =  E5 


(40) 


(41) 

(42) 

(43) 

(44) 


K,  .Q_  +  K,  -Y-  +  K,  _Z„  +  K,  ,Q,  +  K,  +  K,  ,Z_  =  E,  (45) 
6,1  2  6, 22  6, 32  6, 43  6, 53  b,b3  b 

K?  ^2  +  K?  2Y2  +  K?  3Z2  +  K?  4Q3  +  K?  5Y3  +  K?  gZ3  =  E?  (46) 

K8,l92  *  K8,2Y2  +  K8,3Z2  ’■  K8,4Q3  +  K8,5Y3  *  K8,6Z3  =  E8  M7) 


and  for  branch  3 


K20,4Q7  +  K20,5Y7  +  K20,6Z7  =  E20 

K21,1Q7  +  K21,2Y7  +  K21,3Z7  +  K21,4Q8  +  K21,5Y8  +  K21,6Z8  =  E21  (50) 

K22,1Q7  +  K22,2Y7  +  K22,3Z7  +  K22,4Q8  +  K22,5Y8  +  K22,6Z8  =  E22  (51) 

K23,1Q7  +  K23,2Y7  +  K23,3Z7  +  K23,4Q8  +  K23,5Y8  +  K23,6Z8  =  E23  (52) 

K24,1Q8  +  K24,2Y8  +  K24,3Z8  +  K24,4Q9  +  K24,5Y9  +  K24,6Z9  =  E24  (53) 

K25,1Q8  +  K25,2Y8  +  K25,3Z8  +  K25,4Q9  +  K25,5Y9  +  K25,6Z9  =  E25  (54) 

K26,1Q8  +  K26,2Y8  +  K26,3Z8  +  K26,4Q9  +  K26,5Y9  +  K26,6Z9  =  E26  (55) 

Equations  40  and  41  as  well  as  Equations  48  and  49  have  the  form  of  the 
upstream  boundary  Equations  37  and  38,  whereas  Equations  42-47  and 
50-55  are  interior  equations  of  the  form  of  Equations  34-36. 

39.  Equations  40  and  41  with  three  unknowns  can  be  reduced  to 


Q1  =  L1 ,2  +  L1,3Z1 


Y1  =  L2,2  +  L2,3Z1 


in  which  the  coefficient  L's  are  functions  of  K's  and  E's 
stituting  Equations  56  and  57  into  Equations  42-44  then  yields 


E3,3Z1  +  E3,4^2  +  E3,5Y2  +  L3,6Z2  =  M3 


Q_  =  L.  _  +  L,  _Z„ 
v2  4,2  4,3  2 


Y2  =  L5,2  +  L5,3Z2 


The  same  procedure  is  then  used  to  reduce  Equations  45-47  to 


•  •••  ••••••  • 


•  •  •  • 


L6,3Z2  +  L6,4^3  +  L6,5Y3  +  Jj6,6Z3  ”  M6 

(61) 

Q3  =  L7,2  +  L7,3Z3 

(62) 

• 

Y3  =  L8,2  +  L8,3Z3 

(63) 

by  substituting  Equations  58-60  into  Equations  45-47. 

• 

40.  The  same  procedure  is  now  followed  on  branch  3  to 

yield  first 

Q7  =  L19 ,2  +  L19,3Z7 

(64) 

"  • 

Y7  =  L20,2  +  L20,3Z7 

(65) 

then 

• 

L21,3Z7  +  L21,4Q8  +  L21,5Y8  +  L21,6Z8  =  M21 

(66) 

Q8  =  L22,2  +  L22,3Z8 

(67) 

• 

Y8  =  L23,2  +  L23,3Z8 

(68) 

and  finally 

• 

L24,3Z8  +  L24,4Q9  +  L24,5Y9  +  L24,6Z9 

(69) 

Q9  =  L25,2  +  L25,3Z9 

(70) 

• 

Y9  =  L26,2  +  L26,3Z9 

(71) 

41.  Before  proceeding  with  the  forward  sweep  in  branch 

2,  the  junction 

must  be  handled.  The  confluence  Equations  20-23  are  written  as  • 

P28Q4  +  P29Q3  +  P30Q9  =  P1  ^72) 

* 


25 


(73) 


P31Q4  +  P32Y4  +  P33Z4  +  P34Q3  +  P35Y3 


+  P36Z3  +  P37Q9  +  P38Y9  +  P39Z9  ~  P2 


P14Q3  +  P15Y3  +  P16Z3  +  P17Q4  +  P18Y4  +  P19Z4  =  P20 


P21Q9  +  P22Y9  +  P23Z9  +  P24Q4  +  P25Y4  +  P26Z4  =  P27  (75) 

in  which  the  P's  are  functions  of  known  variables.  With  the  Equa¬ 
tions  62,  63,  70,  71,  72,  73,  74,  and  75  a  set  of  eight  linear  equations 
in  nine  unknowns  is  formed.  Two  equations  of  the  form 


Q4  =  L9,2  +  L9,3Z4 


Y4  =  L10,2  +  L10,3Z4 


can  then  be  derived.  The  forward  sweep  continues  on  branch  2  just  as 
performed  previously  on  branches  1  and  3  to  finally  yield 


L14,3Z5  +  L14,4Q6  +  L14,5Y6  *  L14,6Z6  =  M9  (78) 


Q6  =  L15 , 2  +  L15,3Z6 


Y6  =  L16,2  +  L16,3Z6 


42.  Equations  79  and  80  along  with  the  equation  for  the  down¬ 
stream  boundary  of  the  form  of  Equation  39  enables  one  to  compute 
Y6+1  ’  Z6+*  ‘  un^nown  variables  at  the  other  nodes  are  then  com¬ 
puted  by  sweeping  backward,  e.g.  knowing  one  can  compute 

from  Equation  78.  After  computing  ,  Z^  ,  Equation  74  along 

with  Equations  62  and  63  form  a  set  of  three  equations  and  three  unknowns 
from  which  Q3  ,  Y3  ,  Z3  can  be  computed.  Similarly,  Equation  75  along 
with  Equations  70  and  71  enables  the  computation  of  ,  Y^  ,  .  The 

backward  sweep  can  then  be  continued  on  each  branch  until  all  unknown 


variables  have  been  computed.  All  the  coefficients  L's  and  M's  in  the 
equations  can  be  computed  from  recurrence  equations  and  therefore  can  be 
easily  programmed  on  a  digital  computer. 

43.  As  noted  by  Chen,  the  double-sweep  method  offers  the  following 
advantages:  (a)  the  computations  do  not  involve  any  of  the  many  zeros 

in  the  coefficient  matrix  which  saves  considerable  computing  time  and 
(b)  the  computer  storage  required  is  reduced  from  that  required  for  a 
3(IMAX)  x  3(IMAX)  matrix  to  that  of  a  3(IMAX)  x  6  matrix  where  IMAX  is 
the  total  number  of  nodes. 


•  ••♦••••••••••• 
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44.  As  previously  noted,  the  computer  model  is  called  FLOWSED  to 
reflect  the  fact  that  both  flow  and  sediment  movement  are  computed.  The 
basic  equations  solved  in  FLOWSED  are  expressions  of  the  conservation  of 
sediment  mass  and  the  conservation  of  mass  and  momentum  of  the  sediment- 
water  mixture.  Therefore,  the  dependent  variables  computed  are  the  flow 
discharge  and  depth  and  the  bed  elevation,  while  the  independent  vari¬ 
ables  are  time  and  the  distance  along  the  channel.  The  basic  solution 
scheme  is  an  implicit  finite-difference  scheme  in  which  the  resulting 
system  of  linear  algebraic  equations  are  solved  using  the  double-sweep 
algorithm. 


Input  Required 

45.  With  its  present  structure,  all  input  is  read  from  cards  ex¬ 
cept  the  geometry  data  which  can  be  read  from  a  disk  file  if  the  user 
desires.  The  input  data  can  be  broken  into  10  groups;  namely,  general 
control  parameters,  plotting,  branch  information,  junction  information, 
locks  and  dams,  coefficients  in  the  sediment  transport  function  and  eddy 
head  loss  term,  initial  conditions,  information  on  levees,  geometric 
tables,  and  boundary  conditions.  Each  is  briefly  discussed  below.  A 
detailed  list  of  the  input  data  is  presented  in  Appendix  A. 

General  control  information 

46.  Several  parameters  that  describe  the  system  being  modeled  and 
determine  the  various  forms  of  output  must  be  input.  These  include  the 
number  of  net  points,  branches,  junctions,  and  locks  and  dams  as  well  as 
the  numbers  of  the  net  points  at  which  printed  output  is  desired.  Output 
can  be  requested  at  two  different  time-step  intervals  through  the  speci¬ 
fication  of  a  large  and  small  print  interval.  The  small  print  interval, 
of  course,  only  occurs  for  some  specified  length  of  time  during  the 
computations . 

Plotting  data 

47.  CALCOMP  plots  at  as  many  as  20  net  points  can  be  requested. 
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These  plots  can  consist  of  only  a  water-surface  elevation,  discharge,  or 
velocity  time  series  or  various  combinations  of  the  three. 

Branch  information 

48.  With  its  present  dimension  statements,  FLOWSED  can  be  applied 
to  a  system  composed  of  as  many  as  29  branches  containing  a  total  of 
less  than  400  net  points.  For  each  branch,  one  must  input  the  first  and 
last  net  points  on  that  branch  as  well  as  information  related  to  the 
type  of  boundary  condition  to  be  applied  if  the  branch  contains  an  open 
boundary.  It  is  important  to  note  that  all  branches  and  corresponding 
net  points  on  the  main  river  must  be  numbered  before  numbering  the  first 
tributary.  An  example  of  the  proper  numbering  sequence  is  given  below 
for  a  system  composed  of  5  branches  and  29  net  points. 

© 
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Junction  information 

49.  FLOWSED  is  currently  dimensioned  to  handle  14  junctions 
corresponding  to  the  maximum  of  29  branches.  At  each  junction  one  must 
input  the  numbers  of  the  three  branches  composing  that  junction.  Once 
again  it  is  important  to  remember  the  numbering  sequence  of  the  branches. 
Locks  and  dams 

50.  As  many  as  25  locks  and  dams  can  be  handled  by  the  current 
version  of  FLOWSED.  A  descriptive  title  plus  the  pool  elevation  to  be 
maintained  and  the  net  point  immediately  upstream  must  be  input.  In 
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addition,  through  an  input  parameter  one  specifies  whether  either  time- 
varying  water-surface  elevations  or  a  rating  curve  will  be  prescribed 
later. 

Coefficients 

51.  Coefficients  required  in  the  sediment  transport  function  and 
the  eddy  head  loss  term  can  be  prescribed  as  spatially  varying  data  or 
can  be  set  equal  to  constant  values  if  desired.  The  role  of  these  coef¬ 
ficients  can  be  seen  from  the  expressions  presented  for  Se  and  Qg  in 
Equations  10  and  15 . 

Initial  conditions 

52.  Initial  values  of  the  water-surface  elevation,  flow  discharge, 

and  bed  elevation  must  be  input  at  each  net  point.  Simons  et  al.  (1975) 
indicate  that  the  bed  elevation  specified  on  the  upper  Mississippi  River 
application  corresponded  to  the  deepest  1000-ft  width  of  river  channel. 
Initial  values  of  the  discharge  and  depth  are  not  extremely  crucial  since 
their  effect  is  flushed  from  the  system  fairly  quickly.  The  time  re¬ 
quired  for  a  free  surface  wave  traveling  at  a  speed  of  to  traverse 

the  system  twice  is  a  good  rule  of  thumb  for  the  time  required  for 
initial  effects  to  die  out.  If  lateral  inflows  are  specified,  they  must 
also  be  initially  prescribed. 

Levee  information 

53.  One  can  specify  that  certain  reaches  contain  levees  that  will 
be  overtopped.  The  average  elevation  of  the  levee  top  and  the  upstream 
net  point  of  the  levee  must  be  input.  Lateral  outflow  is  then  computed 
from  Equation  24. 

Geometry  data 

54.  The  geometric  tables  constitute  the  majority  of  the  input 
data  required  by  FLOWSED.  At  each  net  point  in  the  system  a  geometric 
table  consisting  of  three  parts  must  be  input.  First,  a  descriptive 
title,  the  river  mileage,  and  the  bank  and  bed  elevations  must  be  input. 
Next,  the  flow  area,  top  width,  and  Manning's  n  versus  elevation  are 
input  for  the  channel.  Finally,  the  floodplain  cross-sectional  area  and 
Manning's  n  versus  elevation  are  input  for  elevations  above  the  top  of 
the  channel.  It  might  be  noted  that  the  river  mileage  of  a  tributary 


must  be  zero  at  its  junction  with  the  main  river.  Mileage  information 
is  used  to  compute  the  computational  spatial  steps,  Ax's  . 

Boundary  conditions 


55.  If  a  branch  has  been  specified  as  one  having  an  open  boundary, 
a  boundary  condition  must  be  prescribed.  At  upstream  boundaries,  the 
flow  and  sediment  discharge  must  be  specified,  although  one  can  request 
that  the  sediment  discharge  be  computed  internally.  At  the  downstream 
boundary  of  the  main  river,  either  a  rating  curve  or  water-surface 
elevations  must  be  input.  Other  required  input  data  are  lateral  inflows 
and  time-varying  pool  elevations.  At  each  time-step,  the  check  on 
time-varying  input  data  follows  the  order  below, 

-  lateral  inflows 

-  all  upstream  boundaries 

-  main  river  downstream  boundary 

-  pool  elevations 


Description  of  Subroutines 


56.  FLOWSED  is  composed  of  a  main  program  that  controls  the  flow 
of  computations  and  16  subroutines  that  perform  various  functions.  A 
brief  description  of  the  role  of  each  is  presented  below. 


LOCKDAM  -  Sets  the  coefficients  in  the  equations  applied 
over  a  computational  reach  containing  a  navigation  lock 
and  dam.  For  example,  to  force  the  flow  discharge  to  be 
the  same  upstream  and  downstream  of  a  dam  the  coefficients 
in  the  equation 


C£,lQi  +  C£,2Yi  +  C£,3Zi  +  C£,4Qi+l  +  C£,5Yi+l  +  C£,6Zi+l  =  E£ 


are  set  to 


£,1  =  1  »  C£t4  =  *  Cl,2  =  C£,3  =  C£,5  =  C£,6  =  E£  =  0  * 
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LINEAR  -  Linearly  interpolates  for  the  flow  area,  top  width, 
and  Manning's  n  from  the  geometric  tables  for  a  particular 
water-surface  elevation.  LINEAR  is  called  many  times  during 
the  computations. 

SEDAREA  -  Computes  the  change  in  the  flow  area  as  the  result  of 
either  erosion  or  deposition. 

FLOOD  -  Determines  the  surface  area  of  the  floodplain  used  in 
computing  q. 

1 

INITIAL  -  Initializes  computations.  Various  minor  computa¬ 
tions  are  performed,  e.g.  the  floodplain  area  which  is  input 
as  a  cross-sectional  area  is  changed  to  surface  area  and 
the  spatial  computational  steps  are  computed. 

JOINTFR  -  Joins  results  from  the  forward  sweep  on  the  main 
river  and  a  tributary  with  the  confluence  equations  so  that 
the  forward  sweep  on  the  main  river  can  then  continue. 

COEFFIT  -  As  noted  in  the  discussion  of  the  solution  scheme, 
the  coefficients  in  the  system  of  linear  algebraic  equations 
are  dependent  upon  information  known  from  the  previous  time 
line.  These  coefficients  are  computed  in  COEFFIT. 

FORWARD  -  Computes  the  coefficients  in  the  forward  sweep  of 
the  double-sweep  solution  algorithm. 

BACKWARD  -  Completes  the  double-sweep  algorithm  by  computing 
the  unknown  variables. 

NEWFLOW  -  Called  at  the  end  of  each  time-step  to  update  flow 
conditions  to  initiate  computations  on  the  next  time  line. 


1 


1 

1 
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DAMRC  -  Computes  pool  elevations  upstream  of  a  lock  and  dam 
from  a  rating  curve  at  the  structure. 

DOWNCOD  -  The  rating  curve  at  the  downstream  end  of  the  main 
river  must  be  cast  into  the  form  of  the  difference  Equa¬ 
tion  39.  This  is  accomplished  by  breaking  the  rating  curve 
up  into  linear  segments.  Each  linear  segment  is  then  de¬ 
fined  by  specifying  the  discharge  and  corresponding  eleva¬ 
tion  at  the  end  of  the  segment  as  well  as  where  the  segment 
crosses  the  elevation  axis.  From  this  information  the 
difference  coefficients  are  computed  in  DOWNCOD. 

BRYCAL  -  Inputs  time -dependent  boundary  data.  Either 
linearly  interpolated  data  or  data  that  have  just  been 
read  are  returned  from  BRYCAL. 

LATERAL  -  Performs  the  same  function  as  BRYCAL  except  that 
lateral  inflows  instead  of  open-boundary  data  are  returned. 

COPOUT  -  Performs  the  function  of  grouping  output  together  as 
well  as  controlling  the  punching  of  restart  cards. 

PLTOUT  -  Generates  CALCOMP  plots,  if  time-history  plots  of 
elevation,  flow  discharge,  or  flow  velocity  are  desired. 


mutation  Cycle 


57.  As  noted,  FLOWSED  is  composed  of  a  main  program  and  16  sub¬ 
routines.  The  main  program  controls  the  computational  cycle  with  the 
various  subroutines  performing  the  tasks  outlined  above.  After  all  input 
data,  except  for  the  time-history  of  lateral  inflows  and  boundary  condi¬ 
tions,  have  been  read  and  INITIAL  has  been  called,  the  basic  loop  in  the 
unsteady  flow  computations  is  entered  in  the  main  program.  Checks  are 
then  made  to  determine  if  and  what  type  of  output  is  to  be  printed. 


Once  the  print  and/or  plot  controls  have  been  exercised,  LATERAL  is 
called  to  return  updated  lateral  inflows  with  a  subsequent  call  to 
BRYCAL  to  return  new  boundary  conditions. 

58.  The  next  step  in  the  computational  cycle  is  either  the  speci¬ 
fication  or  computation  of  the  coefficients  in  the  difference  equations. 
Some  coefficients  are  specified  in  the  main  program  while  those  influ¬ 
enced  by  a  lock  and  dam  are  specified  through  a  call  to  LOCKDAM.  The 
vast  majority  of  the  coefficients  are,  of  course,  computed  from  known 
conditions  through  a  call  to  COEFFIT. 

59.  Once  all  coefficients  have  been  determined,  the  double-sweep 
algorithm  is  initiated.  FORWARD  is  called  for  the  first  branch  on  the 
main  river  and  then  for  the  first  tributary.  JOINTFR  is  then  called  to 
combine  the  results  from  the  two  forward  sweeps  with  the  difference  form 
of  the  confluence  equations.  FORWARD  is  then  called  for  the  main  river 
branch  downstream  of  the  confluence.  This  procedure  is  repeated  until 
the  downstream  end  is  reached.  The  backward  sweep  is  then  initiated  by 
calling  BACKWARD  for  the  last  branch  on  the  main  river  and  then  the 
branch  on  the  main  river  upstream  of  the  last  confluence.  BACKWARD  is 
then  called  for  the  last  tributary.  This  procedure  continues  until  the 
upstream  end  of  the  main  river  is  reached. 

60.  The  final  sequence  of  steps  in  the  computational  cycle  is  the 
calling  of  NEWFLOW  and  a  subsequent  updating  of  the  velocity,  discharge, 
water  depth,  etc.,  arrays  in  the  code.  Control  is  then  transferred  to 
the  beginning  of  the  computational  cycle  and  the  steps  outlined  above 
are  repeated  for  the  next  time-step.  After  all  computations  are  com¬ 
pleted,  COPOUT  and  PLTOUT  are  called  for  the  handling  of  output. 


Output  Furnished 

61.  Both  printed  output  and  time  series  plots  can  be  obtained 
from  FLOWSED.  In  the  printed  output,  either  a  limited  form  of  the 
output  consisting  of  basically  the  water-surface  elevation,  flow  dis¬ 
charge  and  velocity,  bed  elevation,  and  sediment  discharge  at  each  net 
point  requested  is  provided  or  a  more  detailed  listing  of  output  in 


which  much  of  the  geometric  data  are  listed  can  be  obtained.  As  previ¬ 
ously  discussed,  time-history  plots  at  as  many  as  20  net  points  can  be 
obtained.  These  plots  may  consist  of  elevations  and/or  discharges 
and/or  velocities.  Figure  6  is  an  example  of  the  form  of  the  plots  that 
can  be  obtained  from  FLOWSED,  while  an  example  of  printed  output  at  one 
time-step  from  the  limited  output  mode  is  presented  in  Appendix  B. 

Limitations 


62.  FLOWSED  is  a  one-dimensional  unsteady  flow-sediment  implicit 
finite-difference  model  that  can  be  applied  to  a  fairly  general  system 
of  open  channels  containing  navigation  locks  and  dams.  However,  there 
are  limitations  to  its  applicability.  FLOWSED  can  only  be  applied  to  a 
simply  connected  system,  i.e.  closed  loops  within  the  system  cannot  be 
handled.  This  is  illustrated  below. 


Can  be  applied 


Cannot  be  applied 


63.  An  additional  limitation  is  that  the  modeled  system  must 
contain  at  least  one  tributary.  Another  restriction  on  the  physical 
region  is  that  there  can  only  be  one  downstream  boundary.  For  example, 
in  its  present  form,  FLOWSED  cannot  be  directly  applied  to  a  region  such 
as  the  Mississippi  River  near  the  Gulf  of  Mexico  where  flow  enters  the 
Gulf  through  several  distributaries. 

64.  A  major  limitation  of  the  sediment  routing  capability  of 
FLOWSED  is  that  no  provision  for  a  grain-size  distribution  is  allowed. 
The  particular  sediment  transport  function  incorporated  into  FLOWSED  was 
adopted  by  Chen  for  use  on  the  upper  Mississippi  River  and  has  not  been 
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changed.  If  ORD  becomes  more  interested  in  the  sediment  routing  capabil¬ 
ity  of  FLOWSED,  this  is  an  area  that  deserves  additional  investigation. 

65.  In  its  present  form,  there  is  some  limitation  on  the  specifica¬ 
tion  of  boundary  conditions.  At  an  upstream  boundary,  only  flow  dis¬ 
charges  can  be  prescribed;  whereas  at  a  downstream  boundary,  either  a 
rating  curve  or  water-surface  elevations  may  be  specified.  Mathemati¬ 
cally,  one  can  specify  elevations  at  an  upstream  boundary  and/or  dis¬ 
charges  at  a  downstream  boundary  but  some  additional  coding  would  be 
required  to  provide  these  options  in  FLOWSED.  For  ORD's  use  of  the 
model,  these  options  were  not  needed.  It  might  be  noted  that  all  of  the 
limitations  discussed  could  probably  be  removed  without  an  excessive 
effort  being  required  except  for  the  incapability  of  handling  closed 
loops,  which  would  be  a  major  undertaking. 


PART  IV:  DEVELOPMENT  OF  MODEL  INPUT  FOR 
THE  OHIO  RIVER  SYSTEM 


Schematization  of  Physical  Region 

66.  As  noted  in  PART  I,  a  computer  model  called  SOCHMJ  has  been 
applied  to  the  lower  Ohio  River  system  in  previous  studies.  The  first 
application  was  to  the  system  shown  in  Figure  1,  whereas  a  later  applica¬ 
tion  extended  the  model  limits  to  those  shown  in  Figure  2.  The  Ohio 
River  system  to  be  modeled  by  FLOWSED  included  the  previous  area  as  well 
as  the  remainder  of  the  Ohio  River  up  to  Pittsburgh,  Pennsylvania,  and 
portions  of  the  Kentucky,  Great  Miami,  Licking,  Scioto,  Big  Sandy, 
Kanawha,  and  Muskingham  Rivers  as  routed  tributaries.  Figure  7  presents 
a  schematic  of  the  total  system  to  be  modeled,  including  the  location  of 
the  navigation  locks  and  dams.  Table  1  presents  a  description  of  the 
branches  comprising  the  system  while  Table  2  lists  the  pool  elevations 
maintained  by  the  locks  and  dams  in  Figure  7.  It  should  be  noted  that 
the  locks  and  dams  shown  represent  the  system  as  it  existed  in  1981.  In 
the  applications  to  be  presented  later  which  use  historical  data,  some 
of  the  structures  had  not  been  constructed  and  thus  were  not  included  in 
those  applications.  For  example,  the  status  of  the  system  in  May  1974 
is  illustrated  in  Figure  8. 

Development  of  Geometry  Data 


67.  The  geometric  tables  required  by  FLOWSED  were  constructed 
using  data  from  the  Mississippi  Basin  Model  and  cross-sectional  data 
furnished  by  ORD.  The  bulk  of  the  data  downstream  from  McAlpine  Lock 
and  Dam  were  storage  volume  data  collected  from  the  MBM  during  the  early 
1970's.  These  data  were  collected  by  partitioning  the  river  into  small 
reaches  through  the  use  of  sheet-metal  dividers.  Water  was  then  re¬ 
leased  into  the  individual  segments  at  a  known  rate  for  a  known  length 
of  time  with  subsequent  measurements  of  the  elevation  and  water-surface 
top  widths.  Knowing  the  length  of  the  segment,  along  with  an  estimate 
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of  the  top  of  the  bank  elevation,  tables  of  both  channel  and  floodplain 
geometry  could  then  be  constructed. 

68.  While  the  majority  of  the  tables  on  the  lower  Ohio  River  were 
constructed  from  MBM  data,  those  above  McAlpine  Lock  and  Dam  were  con¬ 
structed  from  cross-sectional  data  cards  furnished  by  ORD.  To  convert 
these  data  into  the  tabular  form  required  by  FLOWSED,  a  computer  program 
called  GEOM  was  written.  The  rationale  behind  GEOM  can  best  be  il¬ 


lustrated  by  considering  the  figure  below 


LARGE 

REACH 


Assume  that  cross-sectional  data  have  been  provided  at  each  of  the 
sections  bounding  what  is  referred  to  as  a  small  reach  above.  Further¬ 
more,  assume  that  one  wishes  to  generate  tabular  data  that  are  an  average 
over  a  large  reach.  GEOM  uses  the  trapezoidal  rule  to  compute  the 
channel  area  at  various  elevations  for  each  cross  section.  This  area 
multiplied  by  the  length  of  the  channel  yields  the  channel  volume  in  a 
small  reach  at  a  particular  elevation.  Similarly,  GEOM  also  computes 
the  overbank  areas  and  associated  overbank  volume  at  particular  eleva¬ 
tions  in  the  floodplain  for  each  small  branch.  These  computations  begin 
at  the  downstream  end  and  proceed  up  the  river  until  the  last  small 
reach  is  encountered. 

69.  GEOM  then  checks  to  determine  the  limits  of  each  large  reach. 
Within  each  large  reach,  the  small  reach  volumes  are  summed  at  the 


elevation  entries  associated  with  the  extreme  downstream  section  of  the 
large  reach.  This  of  course  normally  requires  interpolation  among  the 
elevation  entries  of  each  small  reach.  The  total  channel  and  overbank 
volumes  are  then  divided  by  the  length  of  the  large  reach  to  yield  the 
average  cross-sectional  data  of  the  large  reach  at  particular  elevations. 

70.  As  shown  in  Table  1,  12  tributaries  are  modeled  in  the  Ohio 
system.  It  should  be  noted,  however,  that  much  of  the  geometry  data  on 
the  upper  tributaries  was  constructed  from  crudely  estimated  cross 
sections  that  were  often  considered  as  small  rectangular  channels. 

71.  The  complete  Ohio  system  schematized  in  Figure  7  covers  in 
excess  of  1,600  river  miles  and  consists  of  25  branches,  12  junctions, 

20  navigation  locks  and  dams,  and  352  net  points.  These  points  are 
irregularly  spaced  and  are  separated  by  distances  ranging  from  perhaps 
10  miles  on  the  lower  Ohio  and  Mississippi  Rivers  to  5  miles  in  the 
upper  portions  of  the  model.  Of  course,  the  net  points  surrounding 
locks  and  dams  and  tributary-main  river  junctions  are  much  closer 
together. 


Lateral  Inflow  Points 

72.  Lateral  inflows  were  furnished  by  ORD  as  either  gaged  or 
ungaged  data.  The  gaged  flows  represent  flow  from  small  rivers  such  as 
the  Little  Kanawha  which  are  not  handled  as  routed  branches.  The  un¬ 
gaged  data  are  normally  given  as  flow  to  be  distributed  over  relatively 
large  reaches  of  the  Ohio  River,  e.g.  the  Pittsburgh  to  Wheeling  reach. 
Based  upon  an  inspection  of  navigation  maps  showing  natural  drainage 
lines  into  the  Ohio  River,  the  distribution  of  ungaged  lateral  inflows 
presented  in  Table  3  has  been  prescribed.  The  gaged  inflow  points  are 
shown  in  Table  4. 


PART  V:  APPLICATIONS  OF  FLOWSED  TO  THE  OHIO  RIVER  SYSTEM 


73.  FLOWSED  has  been  applied  to  the  complete  Ohio  River  system  in 
a  calibration  mode.  Calibration  runs  were  made  for  historical  floods  in 
1964,  1972,  and  1976.  The  applications  have  been  made  in  a  segmented 
fashion  using  data  furnished  by  ORD  for  different  reaches  of  the  system. 
In  each  of  these  applications,  the  values  for  the  friction  coefficient, 
i.e.  Manning's  n  ,  were  varied  along  the  channel  as  well  as  with  water- 
surface  elevation  until  sufficient  agreement  between  recorded  and  com¬ 
puted  water-surface  elevations  and  flow  discharges  was  realized.  The 
final  n  values  generally  ranged  between  0.025  to  0.035,  and  tended  to 
decrease  with  stage.  In  all  applications,  the  coefficient  C^  in  the 
sediment  transport  function  given  by  Equation  15  was  taken  to  be  1/10  of 
the  value  used  by  Chen  on  the  Mississippi  River  to  minimize  the  effect 
of  the  sediment  computations  on  the  flow.  This  was  done  because  at  the 
present  time  the  primary  interest  in  FLOWSED  is  as  a  flow  model. 

Application  from  Pittsburgh  to  Parkersburg  Using  1972  Data 

74.  Figure  9  shows  that  portion  of  the  Ohio  River  system  extending 
from  Pittsburgh  to  Parkersburg  on  a  larger  scale  than  shown  in  Figure  7. 
As  can  be  seen,  neither  the  Hannibal  nor  the  Willow  Island  Locks  and 
Dams  had  been  constructed  and  thus  the  low-lift  Locks  and  Dams  12-17 
shown  had  to  be  considered.  Table  5  lists  the  pool  elevations  main¬ 
tained  by  these  structures. 

75.  In  analyzing  the  1972  data,  it  was  observed  that  after  con¬ 
trol  was  expected  to  have  been  lost  at  Emsworth  and  Dashields  Locks  and 
Dams,  the  difference  in  the  upstream  and  downstream  elevations  at  each 
structure  indicated  that  the  flow  could  not  be  treated  as  being  free. 

In  discussions  with  personnel  from  ORD  it  was  suggested  that  internal 
rating  curves  be  used  at  these  structures.  These  are  presented  in 
Table  6. 

76.  Figure  10  is  a  plot  of  the  inflow  hydrographs  prescribed  at 
Pittsburgh  and  McConnelsville  along  with  the  lateral  inflow,  while 
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Figure  11  reflects  the  elevation  hydrograph  prescribed  at  Parkersburg. 
Initial  conditions  were  a  steady-state  condition  obtained  by  applying 
FLOWSED  with  the  boundary  conditions  occurring  on  20  June  held  constant. 

77.  Results  from  the  application  are  presented  in  Figures  12-27 
in  which  a  comparison  of  recorded  and  computed  water-surface  elevations 
are  shown.  In  general,  agreement  to  within  1  ft  is  realized.  Again, 
it  should  be  stressed  that  these  results  were  obtained  after  several 
runs  were  made  in  which  Manning's  n  was  varied  from  one  run  to  the 
next.  The  computational  time-step  for  all  applications  was  1  hr. 

Application  from  Parkersburg  to  Huntington 
Using  1964  Data 

78.  Figure  28  is  a  location  map  of  the  Parkersburg  to  Huntington 
reach.  As  can  be  seen  from  a  comparison  of  Figures  7  and  28,  of  the 
current  high-lift  locks  and  dams  only  Gallipolis  existed  in  1964.  In 
addition  to  Gallipolis,  Lock  and  Dam  22  with  a  pool  elevation  of 

551.4'  ft  and  Winfield  Lock  and  Dam  on  the  Kanawha  River  with  an  upstream 
pool  elevation  of  566.0  ft  were  also  in  the  system. 

79.  The  boundary  inflow  hydrographs  at  Parkersburg  and  at 
Charleston  on  the  Kanawha  River  are  presented  in  Figure  29  while  Fig¬ 
ure  30  shows  the  lateral  inflow.  The  lateral  inflow  data  were  varied 
substantially  during  the  calibration  effort  in  order  to  obtain  better 
matching  of  the  recorded  and  computed  water-surface  elevations  or  dis¬ 
charges  presented  in  Figures  31-36.  ORD  personnel  have  indicated  that 
such  a  variation  is  hydrologically  acceptable.  Recorded  elevations  were 
prescribed  as  the  downstream  boundary  condition  at  Huntington.  Agree¬ 
ment  between  computed  and  recorded  elevations  is  generally  in  the  1-ft 
range.  However,  as  illustrated  in  Figure  36,  poor  agreement  was  realized 
in  the  peak  discharge  at  Huntington. 

Application  from  Huntington  to  McAlpine 
Using  1964  Data 

80.  The  Parkersburg  to  McAlpine  reach,  for  which  1964  data  were 
available,  was  broken  into  two  reaches  for  calibration  purposes  with  the 
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belief  that  perhaps  the  calibration  effort  could  be  accomplished  more 
easily  by  treating  each  reach  separately.  Figure  37  is  a  location  map 
for  the  Huntington-McAlpine  reach.  Meldahl  Lock  and  Dam  was  not  in¬ 
cluded  because  the  recorded  elevations  showed  that  the  normal  Meldahl 
pool  elevation  of  485  ft  was  not  maintained  and  thus  the  flow  was  not 
controlled. 

81.  Because  this  application  was  essentially  coupled  with  that  of 
the  Parkersburg  to  Huntington  reach,  computed  flow  discharges  at 
Huntington  were  taken  from  the  previous  application  and  prescribed  as 
the  upstream  boundary  condition  on  the  Ohio.  Recorded  discharges  were 
input  on  the  Big  Sandy,  Scioto,  Licking,  Great  Miami,  and  Kentucky 
Rivers.  The  boundary  inflow  hydrographs  are  presented  in  Figure  38 
while  the  downstream  boundary  elevations  at  McAlpine  are  given  in  Fig¬ 
ure  39.  As  shown  in  the  lateral  inflow  hydrographs  in  Figures  40-43, 
the  ungaged  lateral  inflows  were  varied  as  in  the  previous  application 
in  order  to  increase  the  matching  of  the  computed  and  recorded  water- 
surface  elevations  and  discharges  presented  in  Figures  44-58.  As  in  the 
upstream  applications  discussed  above,  computed  and  recorded  elevations 
generally  agree  to  within  1  ft.  The  one  exception  is  at  the  gage 
upstream  of  Meldahl  Lock  and  Dam.  As  observed  at  Huntington,  poor 
agreement  in  peak  discharge  was  realized  at  Cincinnati  and  at  McAlpine 
Lock  and  Dam. 


Application  from  McAlpine  to  Caruthersville 
Using  1976  Data 

82.  The  lower  Ohio  River  from  McAlpine  Lock  and  Dam  through  its 
junction  with  the  Mississippi  River  to  Caruthersville,  Missouri,  was 
calibrated  using  1976  data.  The  location  map  in  Figure  59  shows  that 
the  Green,  Wabash,  Cumberland,  Tennessee,  and  upper  Mississippi  Rivers 
were  all  handled  as  routed  tributaries.  All  locks  and  dams  that  cur¬ 
rently  exist  except  Smithland  were  in  place  in  1976  and  thus  were  con¬ 
sidered  in  the  application. 

83.  Recorded  flow  discharges  were  input  as  the  boundary  conditions 
at  all  upstream  boundaries  while  the  rating  curve  given  in  Table  7  was 


prescribed  as  the  downstream  boundary  at  Caruthersville.  The  boundary 
inflow  hydrographs  are  presented  in  Figure  60  while  the  lateral  inflow 
hydrographs  are  given  in  Figures  61-65.  Once  again  the  lateral  flows 
were  adjusted  to  achieve  better  matching.  Results  from  the  calibration 
exercise  are  presented  in  Figures  66-78.  Generally,  agreement  between 
recorded  and  computed  results  is  not  as  good  as  in  the  upstream  reaches. 

A  major  reason  is  probably  the  effect  of  the  large  floodplains  that  are 
encountered  in  the  lower  Ohio  River. 

Tributary  Results 

84.  All  results  presented  up  to  now  for  each  of  the  applications 
have  been  only  at  gages  on  the  Ohio  River.  Results  obtained  on  the 
tributaries  were  quite  dependent  on  the  quality  of  the  geometry  data. 

For  tributaries  such  as  the  Big  Sandy,  Cumberland,  Tennessee,  and  upper 
Mississippi  Rivers  on  which  relatively  good  geometry  data  were  available, 
results  were  comparable  to  those  obtained  on  the  Ohio  River,  e.g. 
see  Figures  79-82.  However,  on  tributaries  such  as  the  Scioto  River 
where  the  geometry  data  were  quite  crude,  results  were  not  as  good, 
e.g.  see  Figure  83.  In  addition  to  the  geometry  problem,  average 
daily  values  of  inflow  are  not  sufficient  on  most  of  the  tributaries 
to  properly  represent  the  inflow  hydrograph. 


PART  VI:  SUMMARY  AND  RECOMMENDATIONS 


Summary 

85.  Using  a  computer  code  developed  by  Chen  for  use  on  the  upper 
Mississippi  River  as  the  base,  a  one-dimensional  model  for  open-channel 
unsteady  flow  computations  on  the  complete  Ohio  River  system  has  been 
developed.  As  the  name  FLOWSED  implies,  computations  for  both  flow  and 
sediment  movement  are  made  by  employing  a  sediment  transport  function 
that  is  dependent  upon  the  flow  velocity  and  depth.  FLOWSED  is  a  rela¬ 
tively  general  implicit  finite-difference  model  that  can  be  applied  to  a 
system  containing  any  number  of  tributaries  as  well  as  navigation  locks 
and  dams . 

86.  The  major  limitation  in  FLOWSED 's  use  is  that  it  can  only  be 
applied  to  simply  connected  systems.  In  other  words,  multiple  connected 
systems  containing  closed  loops  cannot  be  handled,  nor  can  distributary 
channels.  Other  limitations  such  as  those  imposed  by  the  particular 
form  of  the  sediment  transport  function  exist  but  could  be  removed 
without  excessive  effort  as  the  need  arises. 

87.  FLOWSED  has  been  applied  to  the  complete  Ohio  River  system, 
broken  into  four  separate  segments,  in  a  calibration  mode.  Historical 
data  from  three  separate  floods  (1964,  1972,  1976)  were  utilized.  The 
results  on  the  Ohio  River  were  excellent  at  some  gages  and  encouraging 
at  others.  However,  results  on  some  of  the  tributaries  where  crude 
geometry  data  were  employed  were  rather  poor. 

Recommendations 

88.  The  normal  procedure  in  the  development  of  a  numerical  model 
for  use  as  a  predictive  tool  is  to  apply  the  model  first  in  a  calibra¬ 
tion  effort  using  data  that  represent  the  system  as  it  currently  exists. 
After  the  model  has  been  calibrated,  a  separate  set  of  data  is  then 
utilized  in  a  verification  application.  This  procedure  has  not  been 
accomplished  in  this  study.  Therefore,  it  is  recommended  that  two  sets 
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of  historical  data  for  the  complete  Ohio  River  system  be  assembled. 

These  data  should  reflect  flows  in  the  system  as  it  currently  exists  if 

possible.  In  addition,  it  should  be  noted  that  based  upon  the  calibra-  • 

tion  effort  to  date,  daily  averaged  inflows  on  the  tributaries  as  well 

as  at  the  upper  limit  of  the  Ohio  River  at  Pittsburgh  are  inadequate 

(i.e.,  flows  should  be  given  a  shorter  time  interval,  preferably  hourly). 

89.  After  assembling  the  two  sets  of  historical  data,  FLOWSED  • 

should  be  recalibrated  using  one  set  of  data  and  verified  using  the 

second  set.  If  refined  geometry  data  are  available  on  the  tributaries, 
these  should  be  employed.  In  this  effort  it  is  suggested  that  FLOWSED 
be  applied  to  the  complete  system  rather  than  in  a  segmented  fashion  as  • 

was  done  in  the  present  study. 

90.  Other  recommendations  center  around  suggested  improvements  to 
the  computer  code.  To  better  handle  flows  through  river  bends  and 

floodways,  e.g.  the  New  Madrid  Floodway  near  the  Ohio-Mississippi  junc-  • 

tion,  work  done  at  SOGREAH  in  France  (1978)  should  be  investigated. 

Such  work  centers  around  the  use  of  two-dimensional  cells  in  floodplains 
which  are  connected  with  the  normal  one-dimensional  channel  computations. 

In  addition,  minor  limitations  of  the  current  version  of  the  model  such  • 

as  forcing  the  modeled  system  to  contain  at  least  one  tributary  should 
be  removed. 

• 


REFERENCES 


Chen,  Y.  H.  1973  (Mar).  "Mathematical  Modeling  of  Water  and  Sediment 
Routing  in  Natural  Channels,"  Ph.  D.  Dissertation  at  Colorado  State 
University,  Fort  Collins,  Colo. 

Fread,  D.  L.  "Theoretical  Development  of  Implicit  Dynamic  Routing 
Model,"  presented  at  Dynamic  Routing  Seminar,  Lower  Mississippi  River 
Forecast  Center,  Slidell,  La.,  13-17  Dec  1976. 

Garrison,  M.  J. ,  Granju,  J.-P.,  and  Price,  T.  J.  1969  (Sep).  "Unsteady 
Flow  Simulation  in  Rivers  and  Reservoirs — Applications  and  Limitations," 
Journal,  Hydraulics  Division,  American  Society  of  Civil  Engineers, 


Vol  95,  No.  HY5,  p  1559;  presented  at  ASCE  Hydraulics  Division  Specialty 
Conference;  Cambridge,  Mass.,  21-23  Aug  1968. 

Johnson,  B.  H.  1974  (Sep).  "Unsteady  Flow  Computations  on  the  Ohio- 
Cumber! and-Tennessee-Mississippi  River  System,"  Technical  Report  H-74-8, 
U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE,  Vicksburg,  Miss. 

_ .  1977  (Oct).  "A  Mathematical  Model  for  Unsteady-Flow  Compu¬ 
tations  Through  the  Complete  Spectrum  of  Flows  on  the  Lower  Ohio  River," 
TR-77-18,  U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE, 
Vicksburg,  Miss. 

Johnson,  B.  H.,  and  Senter,  P.  K.  1973  (Jun).  "Flood  Routing  Procedure 
for  the  Lower  Ohio  River,"  MP-H-73-3,  U.  S.  Army  Engineer  Waterways 
Experiment  Station,  CE,  Vicksburg,  Miss. 

Simons,  D.  B. ,  et  al.  1975  (Jul).  "Environmental  Inventory  and  Assess¬ 
ment  of  Navigation  Pools  24,  25,  and  26,  Upper  Mississippi  and  Lower 
Illinois  Rivers,"  CR-Y-75-3,  U.  S.  Army  Engineer  Waterways  Experiment 
Station,  CE,  Vicksburg,  Miss. 

SOGREAH  Consulting  Engineers.  1978  (Oct).  "Modelling  of  Unsteady  Flow 
in  River  and  Flood  Plain  Networks  Using  the  CARIMA  System,"  In-house 
report,  France. 


Table  1 

Schematization  of  the  Ohio  System 


Branch 

No. 

First 

Net  Point 

Last 

Net  Point 

Description 

1 

1 

63 

Ohio  R.  from  Pittsburgh  to  Muskingham  R. 

2 

64 

85 

Ohio  R.  from  Muskingham  R.  to  Kanawha  R. 

3 

86 

100 

Ohio  R.  from  Kanawha  R.  to  Big  Sandy  R. 

4 

101 

111 

Ohio  R.  from  Big  Sandy  R.  to  Scioto  R. 

5 

112 

140 

Ohio  R.  from  Scioto  R.  to  Licking  R. 

6 

141 

145 

Ohio  R.  from  Licking  R.  to  Great  Miami  R. 

7 

146 

157 

Ohio  R.  from  Great  Miami  R.  to  Kentucky  R. 

8 

158 

194 

Ohio  R.  from  Kentucky  R.  to  Green  R. 

9 

195 

205 

Ohio  R.  from  Green  R.  to  Wabash  R. 

10 

206 

217 

Ohio  R.  from  Wabash  R.  to  Cumberland  R. 

11 

218 

220 

Ohio  R.  from  Cumberland  R.  to  Tennessee  R. 

12 

221 

230 

Ohio  R.  from  Tennessee  R.  to  Mississippi  R 

13 

231 

243 

Lower  Mississippi  R.  to  Caruthersville ,  Mo 

14 

244 

257 

Muskingham  R. 

15 

258 

271 

Kanawha  R. 

16 

272 

280 

Big  Sandy  R. 

17 

281 

295 

Scioto  R. 

18 

296 

305 

Licking  R. 

19 

306 

315 

Great  Miami  R. 

20 

316 

324 

Kentucky  R. 

21 

325 

331 

Green  R. 

22 

332 

338 

Wabash  R. 

23 

339 

342 

Cumberland  R. 

24 

343 

346 

Tennessee  R. 

25 

347 

352 

Upper  Mississippi  R. 

Table  2 

Locks  and  Dams  in  Ohio  System 


Lock  and  Dam 

Pool  Elevation 
ft  NGVD 

Upstream 
Net  Point 

Emsworth 

710.0 

4 

Dashields 

692.0 

9 

Montgomery 

682.0 

17 

New  Cumberland 

664.5 

25 

Pike  Island 

644.0 

33 

Hannibal 

623.0 

46 

Willow  Island 

602.0 

58 

Belleville 

582.0 

70 

Racine 

560.0 

78 

Gallipolis 

538.0 

89 

Greenup 

515.0 

108 

Meldahl 

485.0 

133 

Markland 

455.0 

153 

McAlpine 

420.0 

167 

Cannelton 

383.0 

182 

Newburgh 

358.0 

190 

Uniontown 

342.0 

202 

Smithland 

320.0 

214 

Lock  &  Dam  52 

302.0 

225 

Winfield 

566.0 

263 

i 


Table  3 

Ungaged  Lateral  Inflow  Distribution 


Net  Percent 


Point 

Ungaged  Reach 

of  Flow 

3 

Pittsburgh-Wheeling 

5 

8 

Pittsburgh-Wheeling 

5 

10 

Pittsburgh-Wheeling 

5 

13 

Pittsburgh-Wheeling 

5 

16 

Pittsburgh-Wheeling 

7 

20 

Pittsburgh-Wheeling 

13 

24 

Pittsburgh-Wheeling 

13 

28 

Pittsburgh-Wheeling 

17 

32 

Pittsburgh-Wheeling 

30 

100 

36 

Wheeling-St.  Mary 

18 

40 

Wheeling-St.  Mary 

13 

45 

Wheeling-St.  Mary 

16 

51 

Wheeling-St.  Mary 

30 

55 

Wheeling-St.  Mary 

23 

100 

67 

St.  Mary-Pomeroy 

10 

68 

St.  Mary-Pomeroy 

10 

73 

St.  Mary-Pomeroy 

30 

76 

St.  Mary-Pomeroy 

20 

77 

St.  Mary-Pomeroy 

10 

79 

St.  Mary-Pomeroy 

20 

100 

81 

Pomeroy-Huntington 

20 

83 

Pomeroy-Huntington 

10 

84 

Pomeroy-Huntington 

10 

87 

Pomeroy-Huntington 

10 

91 

Pomeroy-Huntington 

20 

97 

Pomeroy-Huntington 

30 

100 

105 

Huntington-Maysville 

30 

114 

Huntington-Maysville 

10 

120 

Huntington-Maysville 

30 

124 

Huntington-Maysville 

30 

Net  Percent 

Point  Ungaged  Reach  of  Flow 

129  Maysville-Cincinnati  30 

132  Maysville-Cincinnati  40 

138  Maysville-Cincinnati  30 

100 

147  Cincinnati-Markland  20 

151  Cincinnati-Markland  80 

100 

160  Markland-McAlpine  50 

165  Markland-McAlpine  50 

100 

174  McAlpine-Evansville  22 

177  McAlpine-Evansville  12 

180  McAlpine-Evansville  33 

189  McAlpine-Evansville  33 

100 

201  Evansville-Golconda  50 

210  Evansville-Golconda  50 

100 

213  Golconda-Metropolis  14 

215  Golconda-Metropolis  14 

219  Golconda-Metropolis  14 

223  Golconda-Metropolis  58 

100 

225  Metropolis-Cairo  50 

228  Metropolis-Cairo  50 

100 

347  Thebes-Cairo  50 

349  Thebes-Cairo  50 


100 


'  • 


100 


Table  4 

Gaged  Inflow  Points 


Net  Point 

Description  of  Inflow 

13 

Beaver  River 

66 

Little  Kanawha  River 

69 

Hocking  River 

96 

Guyandot  River 

99 

Twelve  Pole  Creek 

106 

Little  Sandy  River 

111 

Little  Scioto  River 

140 

Little  Miami  River 

171 

Salt  River 

• 

Table  5 

Pool  Elevations  Maintained  by  Low-Lift 
Structures  in  Pittsburgh-Parkersburg  Reach  • 

Pool  Elevation  Maintained 
Structure  ft  NGVD 

Lock  &  Dam  12 

Lock  &  Dam  13 

Lock  &  Dam  14 

Lock  &  Dam  15 

Lock  &  Dam  16 

Lock  &  Dam  17 


626.5 

616.2 

610.7 

602.3 

594.3 
586.9 


Table  6 

Rating  Curves  Applied  at  Emsworth  and  Dashields 


1 


EMSWORTH 


DASHIELDS 


Elevation 

Discharge 

Elevation 

Discharge 

ft  NGVD 

cfs 

ft  NGVD 

cfs 

707.7 

37,500 

695.25 

32,500 

709.2 

56,000 

698.1 

110,000 

711.3 

175,000 

700.5 

170,000 

• 

721.5 

385,000 

713.5 

385,000 

715.0 

280,000 

708.0 

300,000 

713.0 

222,000 

703.75 

215,000 

711.6 

180,000 

700.3 

160,000 

• 

710.0 

152,500 

699.3 

140,000 

709.9 

130,000 

695.25 

32,500 

709.0 

117,500 

708.0 

100,000 

• 

707.75 

75,000 

Table  7 

Rating  Curve  Employed  at  Caruthersville* 


Elevation 

Discharge 

ft  NGVD 

cfs 

255.0 

200,000 

• 

257.0 

320,000 

259.0 

440,000 

261.00 

575,000 

263.0 

700,000 

• 

265.0 

820,000 

267.0 

940,000 

269.0 

1,050,000 

271.0 

1,170,000 

• 

273.0 

1,290,000 

275.0 

1,410,000 

277.0 

1,550,000 

279.0 

1,710,000 

• 

281.0 

1,920,000 

283.0 

2,240,000 

285.0 

2,900,000 

• 

*  Obtained  from  data 
U.  S.  Army  Engineer 

provided  by  the 

District,  Memphis. 

•  ••••••••••• 
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Figure  2.  Physical  limits  of  second  application  of  SOCHMJ 
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Figure  3.  Ohio  River  Basin 
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Figure  6.  Typical  example  of  FLOWSED's  plotting  capability 
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Figure  7.  Location  map  for  tributaries,  and  locks 
and  dams  on  the  Ohio  River  model 


Status  of  lock  and  dam  construction  on  the  Ohio  River,  May  1974 
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Figure  9.  Location  map  for  Pittsburgh-Parkersburg  reach 
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Figure  10.  Inflow  hydrographs  for  Pittsburgh  to  Parkersburg  application 


Figure  11.  Boundary  elevations  prescribed  at  Parkersburg  in  Pittsburgh-Parkersburg  application 
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Figure  12.  Recorded  versus  computed  elevations  at  Pittsburgh 


Figure  13.  Recorded  versus  computed  elevations  at  Emsworth  L&D,  upper  gage 


722 


Recorded  versus  computed  elevations  at  Dashields  L&D,  upper  gage 


Recorded  versus  computed  elevations  at  Montgomery  L&D,  lower  gage 


»*2 


Figure  19.  Recorded  versus  computed  elevations  at  New  Cumberland  L&D,  upper  gage 


•72 


Figure  20.  Recorded  versus  computed  elevations  at  New  Cumberland  L&D,  lower  gage 


654 


b 


•  •  •  •  • 


Figure  22.  Recorded  versus  computed  elevations  at  Pike  Island  L&D,  lower  gage 


652 


Figure  23.  Recorded  versus  computed  elevations  at  L&D  12,  upper  gage 
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Figure  25.  Recorded  versus  computed  elevations  at  L&D  15,  upper  gage 


Figure  26.  Recorded  versus  computed  elevations  at  L&D  16,  upper  gage 


'arkersburg-Huntington  application 


Figure  30.  Ungaged  lateral  inflow  from  Parkersburg  to  Huntington 
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Recorded  versus  computed  elevations  at  Parkersburg 


Figure  32.  Recorded  versus  computed  elevations  at  Dam  22,  upper  gage 
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Recorded  versus  computed  elevations  at  Point  Pleasant 


Recorded  versus  computed  elevations  at  Gallipolis  L&D,  upper  gage 
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Recorded  versus  computed  discharge  at  Huntington 
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Figure  37.  Location  map  for  Huntington-VcAlpine  reach 


Figure  38.  Inflow  boundary  conditions  for  Huntington  to  McAlpine  application 


*50 
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Figure  40.  Ungaged  lateral  inflow  in  Huntington-Maysville  reach 
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Figure  44.  Recorded  versus  computed  elevations  at  Huntington 


Recorded  versus  computed  elevations  at  Greenup  L&D,  lower  gage 


COMPUTED 


Recorded  versus  computed  elevations  at  Portsmouth 


Recorded  versus  computed  elevations  at  Maysville 
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Recorded  versus  computed  elevations  at  Meldahl  L&D,  upper  gage 


Figure  52.  Recorded  versus  computed  elevations  at  Meldahl  L&D,  lower  gage 


Recorded  versus  computed  elevations  at  Cincinnati 
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Figure  56.  Recorded  versus  computed  elevations  at  Markland  L&D,  upper  gage 


475 


Figure  57.  Recorded  versus  computed  elevations  at  Markland  L&D,  lower 


Recorded  versus  computed  discharges  at  McAlpine  L&D 
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Figure  59.  Location  map  for  McAlpine  to  Caruthersville  application 


Figure  60.  Boundary  inflow  hydrograph  for  McAlpine-Caruthersville  application 


Figure  61.  Ungaged  lateral  inflow  in  McAlpine-Evansville  reach 


Figure  62.  Ungaged  lateral  inflow  in  Evansville-Golconda  reach 


Figure  66.  Recorded  versus  computed  elevations  at  McAlpine  L&D 


Figure  69.  Recorded  versus  computed  elevations  at  Newburgh  L&D,  upper  gage 


Figure  70.  Recorded  versus  computed  elevations  at  Newburgh  L&D,  lower 


Figure  71.  Recorded  versus  computed  elevations  at  Uniontown  L&D 


Figure  72.  Recorded  versus  computed  elevations  at  Uniontown  L&D,  lover  gage 


Figure  73.  Recorded  versus  computed  elevations  at  L&D  51 


Recorded  versus  computed  discharge  L&D  51 


Figure  75.  Recorded  versus  computed  elevations  at  Paducah 
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Figure  76.  Recorded  versus  computed  elevations  at  L&D  52 


Recorded  versus  computed  discharge  at  L&D  52 


versus  computed  elevations  at  Cairo 
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Figure  79.  Recorded  versus  computed  elevations  at  Louisa  (Big  Sandy  River) 


Figure  80.  Recorded  versus  computed  elevations  at  Barkley  L&D  (Cumberland  River) 


Recorded  versus  computed  elevations  at  Thebes  (upper  Mississippi  River) 


sus  computed  elevations  at  Higby  (Scioto  River 


APPENDIX  A:  LIST  OF  FLOWSED  INPUT  CARDS 


1.  TITLE (I) 


(10A8) 


2.  IGEOM,  ILUG,  I  PUNCH,  I BACK 


(415) 


IGEOM  =  0  -  Geometry  tables  printed 
IGEOM  =  1  -  Tables  not  printed 

ILUG  *  Unit  from  which  geometry  tables  will  be  read 
IPUNCH  *  1  -  Restart  cards  will  be  punched 
IPUNCH  =  0  -  No  restart  cards 


I  BACK 


0  -  Less  printed  output 


3.  NC,  NBRS ,  NJUNC ,  NDAMS ,  NXMAIN,  ISTAGE,  IFPLN,  TOTALT,  TSTEP  (715, 

10X,  2F10.0) 


NJUNC 


=  Total  number  of  net  points  minus  one 
■  Total  number  of  branches 
-  Total  number  of  junctions  (-1) 


NDAMS  *  Total  number  of  dams 


NXMAIN  = 
ISTAGE  - 
IFPLN  = 
TOTALT  * 
TSTEP  = 


Last  net  point  on  main  river 

Number  of  entries  in  channel  geometry  tables 

Number  of  entries  in  floodplain  tables 

Number  of  days  of  computations 

Time-step  in  seconds 


4.  MNTH,  KDAY,  KYEAR,  TIME,  NOXS 


(315,  5X,  F10.0,  15) 


MNTH  -  Starting  month 

KDAY  -  Starting  day 

KYEAR  *  Starting  year 

TIME  *  Starting  time  on  a  24-hr  clock 

NOXS  ■  Number  of  stations  at  which  plots  are  desired 


9 


5.  NS TAT,  IPRINT,  INTVG,  INTVD,  INTVP  (515) 

NSTAT  =  Number  of  net  points  at  which  output  is  desired 
IPRINT  *  0  -  Limited  output 
IPRINT  *  1  -  Detailed  output 
INTVG  *  Major  print  interval 

INTVD  ■  Print  interval  for  particular  days  (see  next  card) 

INTVP  ■  Interval  for  placing  points  in  plot  file 
(If  INTVD  ■  INTVG,  insert  blank  card  for  card  6) 

6.  STOP,  ISDDP,  ISMDP,  ISYDP,  ETDP ,  IEDDP ,  IEMDP ,  IEYDP  (F10.0,  315,  5X, 

F10.0,  315) 

STOP  ■  Starting  time  on  24-hr  clock  for  small  print  interval 

ISDDP  *  Starting  day  for  small  print  interval 

ISMDP  -  Starting  month  for  small  print  interval 

ISYDP  ■  Starting  year  for  small  print  interval 

ETDP  -  Ending  time  on  24-hr  clock  for  small  print  interval 

IEDDP  «  Ending  day  for  small  print  interval 

IEMDP  *  Ending  month  for  small  print  interval 

IEYDP  -  Ending  year  for  small  print  interval 

7.  (NPRINT(I) ,  1-1,  NSTAT)  (1615) 

NPRINT  ■  Net  point  numbers  at  which  output  is  requested 

8.  IPLT,  ISCL,  IYCN,  I0P1  (415) 

IPLT  -  0  -  No  plots 

IPLT  •  1  -  Elevation  plots 

IPLT  -  2  -  Discharge  plots 

IPLT  ■  3  -  Elevation,  discharge  plots 

IPLT  *  4  -  Velocity  plots 

IPLT  -  5  -  Elevation,  velocity  plots 


A2 


IPLT  ■  6  -  Discharge,  velocity  plots 

IPLT  *  7  -  Elevation,  discharge,  velocity  plots 

ISCL  *  Number  of  minutes/ inch  on  x-axis 

IYLN  *  Length  of  y-axis  in  inches  (0  defaults  to  8  in.) 

I0P1  ■  0  -  Different  y-interval  on  each  plot 

I0P1  ■  1  -  Same  y-interval  on  each  plot 

9.  (STITLE(II,I) ,  II-l,  4),  NPLOT(l)  (4A8,  3X,  15) 

STITLE  -  Title  of  station 

NPLOT  ■  Net  point  number 

(Repeat  card  9  NOXS  times) 

10.  ID,  IBRNCH(I.l),  IBRNCH(I,2) ,  IBRS(I),  ISED(I) ,  IBC(I)  (615) 

ID  ■  Branch  number 

*IBRNCH(I,1)  *  First  net  point  on  branch 
IBRNCH(I ,2)  *  Last  net  point  on  branch 
IBRS(I)  -  1  -  Branch  has  upstream  outer  boundary 
IBRS(I)  ■  0  -  Branch  is  an  interior  branch 
IBRS(I)  -  -1  -  Branch  has  downstream  outer  boundary 
ISED(I)  *  0  -  Sediment  at  boundary  computed  internally 

ISED(I)  *  1  -  Sediment  at  boundary  input 

IBC(I)  *  -1  -  Rating  curve  will  be  used  if  this  branch  con¬ 

tains  downstream  outer  boundary 

IBC(I)  *  0  -  If  this  is  an  interior  branch  or  discharge  is  Input 

at  its  boundary 

IBC(I)  *  1  -  Elevations  will  be  input  if  this  branch  con¬ 

tains  downstream  outer  boundary 

(There  will  be  one  card  10  for  each  branch.  All  mainstem  branches  should 
be  numbered  before  numbering  tributaries.) 


11.  ID,  (IJUN(I,K) ,K“1,3) ,  AL(J) ,  AL(J+1) 
ID  *  Junction  number 


(415,  10X,  2F10.0) 


IJUN(I,1)  ■  Number  of  upstream  branch  on  main  river 

IJUN(I,2)  -  Number  of  tributary  branch 

IJUN(I,3)  »  Number  of  downstream  branch  on  main  river 

AL(J)  »  Velocity  head  correction  factor  associated  with  junc¬ 
tion  of  upstream  main  river  and  downstream  main  river 

AL(J+1)  -  Velocity  head  correction  factor  associated  with  junc¬ 

tion  of  tributary  and  downstream  river 

(There  will  be  one  card  11  for  each  junction.) 


12.  TDAM(I),  HSET(I),  NL(I),  NVARY(l) 


(A8,  2X,  F10.0,  215) 


HSET(I) 

NL(I) 

NVARY(I) 

NVARY(l) 

NVARY(I) 

13.  ISCON 


Description  of  dam 
Elevation  maintained  by  dam 
Net  point  immediately  upstream  of  dam 
0  -  Normal  dam 

1  -  Time-varying  elevations  of  pool  will  be  input 

2  -  Rating  curve  will  be  input  for  this  dam 


ISCON  -  0  -  Sediment  coefficients  will  be  read  in  as  constants 

ISCON  -  1  -  Sediment  coefficients  will  be  spatially  varying 

(If  ISCON  -  1,  skip  card  14) 


14.  CRT,  CMS,  CNS 


(3F10.0) 


See  definitions  below. 


(If  ISCON 


0,  go  to  card  18) 


15.  (CKT(I) ,1*1, NX) 


(8F10.0) 


(NX  *  Total  number  of  net  points) 

CKT(I)  -  Sediment  transport  coefficient  in  C  -  kvcmdcn  where 
CKT  ■  K  .  A  value  between  0.000005  ind  0.000015  is  ap¬ 
propriate  in  upper  Mississippi  River. 


A4 


1 


16.  ( CMS ( I ) , 1=1 , NX) 


(8F1O.0) 


CMS ( I ) 


=  CM  in  the  equation  above.  Value  of  3.4  is  appropriate 
upper  Mississippi  River. 


17.  (CNS(I) ,  1=1,  NX) 


(8F10.0) 


CNS(I) 


=  CN  in  the  equation  above.  Value  of  -1.0  is  appropriate 
in  upper  Mississippi  River. 


18.  CP 


(F10.0) 


19.  IEDYHD 


=  Initial  volume  of  sediment  per  unit  volume  of  bed  ma¬ 
terial 


IEDYHD 
IEDYHD 
(If  IEDYHD 


=  0  -  Eddy  head  loss  coefficients  set  to  0 
=  1  -  Eddy  head  loss  coefficients  read  in 
=  0,  Skip  card  20) 


20.  (CKE(I),  1=1, NX) 

(NX  =  Total  number  of  net  points) 

CKE(I)  =  Coefficient  in  eddy  loss  term 


(8F10.0) 


21.  (H0(I) ,1=1 ,NX) 


(8F10.0) 


H0(I)  =  Initial  water-surface  elevation  in  ft  at  each  net  point 


22.  (Q ( I ) , 1) , 1=1 ,NX) 


(8F10.0) 


Q(I,1)  =  Initial  discharge  in  cfs  at  each  net  point 


23.  NRCH, (IRCH(I) , 1=1 ,NRCH) 


(1615) 


! 


I 

8 

1 

4 

'  4 


4 


4 


NRCH 
IRCH(I) 
(If  NRCH 


■  Total  number  of  reaches  containing  lateral  Inflows 
*  Numbers  of  the  reaches  containing  lateral  inflow 

■  0,  Skip  the  next  card) 


24.  C0NCL,(QL2P(K),K-1,  NRCH) 


(8F10.0) 


CONCL 

QL2P(K) 


3  3 

*  Concentration  of  sediment  in  ft  /ft  contained  in  the 
lateral  inflow 

*  Lateral  inflow  at  the  initial  time  line  in  cfs 


25.  NLEVEE,  (ILEVEE(I) ,I«1, NLEVEE) 


(1615) 


NLEVEE  «•  Number  of  reaches  with  levees 
ILEVEE(I)  *  Upstream  net  points  of  reaches  with  levees 
(If  NLEVEE  *  0,  Skip  next  card) 


26.  (ELEVEE(K),  K-l, NLEVEE)  (8F10.0) 

ELEVEE(K)  *  Average  elevation  of  the  top  of  the  levee  along  this 
reach  in  feet 


27.  RANGE (I) ,  XL(I),  DUM,  ZF(I) ,  Z0(I) ,  BETA(I)  (A8, 2X, 5F10.0) 


RANGE(I) 

XL(I) 

DUM 

ZF(I) 

ZO(I) 

BETA(I) 


*  Description  of  I'th  net  point 

*  River  mileage  of  I'th  net  point.  Tributary  mileage  is 
zero  at  the  junction 

■  Space  for  top  width.  Can  leave  blank 

*  Top  bank  elevation  for  I'th  net  point 

*  Bed  elevation  of  I'th  net  point 

*  Momentum  correction  factor 


28.  HI(I , J) ,  AI (I , J) ,  TI(I,J),  RNI(I,J) 


(4F10.0) 


HI  ■  Elevation  of  channel  geometry 

AI  ■  Flow  area 

TI  ■  Top  width 


•  •  • 


•  •  • 


A6 


's 

• 

RNI 

-  Manning ' s  n 

(Card  28  is  repeated  ISTAGE  times) 

( 

29.  HF(I,J) , 

AFI(I,J),  RNIFP(I,J)  (3F10.0) 

• 

: 

HF 

*  Elevation  of  floodplain  geometry 

AFI 

*  Cross-sectional  area  of  the  floodplain  at  elevation 

HF(I, J) 

• 

RNIFP 

*  Manning  ' s  n 

* 

(Card  29  is  repeated  IFPLN  times) 

ni 

(Cards  27,  28,  and  29  are  repeated  in  sequence  for  each  net  point.  They 
are  read  from  unit  ILUG) 

-  •  “ 

(If  IBC(I) 

*  1,  Cards  30  and  31  are  omitted) 

> 

\ 

30.  NS  EG 

(15) 

• 

NS  EG 

*  Number  of  linear  segments  approximating  the  rating  curve 

# 

31.  QRC(I) , 

HRC (I) ,  SEE(I)  (3F10.0) 

* 

QRC 

*  Discharge  at  the  end  of  the  I'th  linear  segment 

•  ' 

HRC 

»  Elevation  of  water  surface  corresponding  to  the  end  of  the 

I'th  linear  segment 

SEE 

*  Elevation  of  channel  bottom  -  intercept  of  this  segment  on 
the  elevation  axis. 

m 

l 

(Card  31  is 

repeated  for  each  linear  segment) 

(If  NVARY(I) 

+  2,  Omit  card  32) 

.  i 

.  • 

32.  KRC(I) ,  QLIMIT(I) ,  QCHECKO(I),  ODRCF(I), 

(QDRC(K,I),HDRC(K,I),K-1,KRC(I))  (I5.7F10.0/8F10.0) 

• 

.  i 

KRC 

■  Number  of  entries  in  rating  curve  table  at  the  I'th  dam 

- 

QLIMIT 

■  Discharge  below  which  a  fixed  water-surface  elevation  is 
prescribed 

• 

r- 

QCHECKO 

*  Initial  discharge  of  I'th  dam 

QDRCF 

*  Discharge  above  which  the  falling  portion  of  the  rating 
curve  will  be  used  if  the  discharge  is  decreasing 

; 

A7 

• 

* 

•  •  '  • 

•  •••••••••• 

•  • 

QDRC  -  Discharge  In  rating  curve  table 
HDRC  ■  Water-surface  elevation  corresponding  to  QDRC 
(Card  32  Is  repeated  for  each  dam  vlth  NVARY  -  2) 

(Card  33  is  omitted  If  NRCH  -  0) 

33.  J,  (QL2(K) .K-l.NRCH)  (I5,5X,7F10.0/8F10.0) 

J  -  Number  of  time-steps  before  new  lateral  Inflows  will  be 

Input 

QL2  -  Lateral  inflow  In  cfs 

34.  Q(I ,2) ,  IQCK(I)  (F10.0.I5) 

Q  -  New  boundary  discharge  in  cfs 

IQCK  -  Number  of  time-steps  before  a  new  discharge  will  be  input 

35.  QS(1,2),  IQSCK(I)  (F10.0.I5) 

QS  ■  New  sediment  discharge  at  boundary  In  cfs 

IQSCK  -  Number  of  time-steps  before  a  new  sediment  discharge  will 
be  Input 

(Cards  34  and  35  are  repeated  for  each  upstream  boundary) 

(If  IBC(I)  4  1,  Omit  card  36) 

36.  H(I,2),  IHCK(I)  (F10.0,I5) 

H  “  New  downstream  boundary  water  surface  elevation 

IHCK  -  Number  of  time-steps  before  a  new  elevation  will  be  Input 

(If  NVARY(K)  -  1  for  K  -  1,  NDAMS  omit  card  37) 

37.  HSET(K) ,  IHSET(K)  (FIO.0,15) 

HSET  -  New  water-surface  elevation  to  be  set  upstream  of  the  K'th 
dam 

«  Number  of  time-steps  before  a  new  HSET  Is  read 


IHSET 


(Repeat  card  37  for  each  dam  with  NVARY  =  1) 

Remember,  at  each  time-step  the  check  is  on  lateral  inflows  first,  then 
on  boundary  conditions,  and  finally  on  time-varying  pools. 


•  '• 
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APPENDIX  C:  NOTATION 


< 


A 


C1’C2,C3 

E 

f 

g 

H 

h 

K 

K 

e 

L 

M 

n 

P 

P 

q 


Total  cross-sectional  area  of  channel 
Cross-sectional  area  of  floodplain 

Derivative  of  A  with  respect  to  channel  distance  at  a 
constant  flow  depth 

Top  width  of  water  surface 

Average  suspended  sediment  concentration 

Wash-load  concentration 

Coefficients  in  the  sediment  transport  function 
Coefficient  in  difference  equations 
General  function 

Acceleration  due  to  gravity,  general  function 

Height  of  water  over  a  levee 

Water-surface  elevation 

Coefficient  in  difference  equations 

Coefficient  in  eddy  head  loss  term 

Coefficient  in  solution  scheme 

Coefficient  in  solution  scheme 

Coefficient  in  Manning's  equation 

Porosity  of  bed  layer 

Wetted  perimeter  of  a  section 

Unit  lateral  outflow  from  levee  overtopping 

Unit  lateral  inflow  rate  of  water 

Unit  lateral  outflow  rate  of  water  to  the  floodplain 

Unit  lateral  inflow  rate  of  water  from  a  tributary 

Unit  lateral  inflow  rate  of  sediment 

Unit  lateral  inflow  rate  of  sediment  from  a  tributary 

Unit  lateral  outflow  rate  of  sediment  to  the  floodplain 

Flow  discharge 
Sediment  discharge 
Wash-load  discharge 


Cl 


* 


*  ' 


•  ' 


•  • 


R 


1 


1 

\ 


H 

i- 

i. 


! 


S 

e 


t 

V 


x 


y 

Y 

z 

Ah 
At 
A Z 
AZf 

a 

P 

P 

Ps 

Pw 

To 

a/at 

a/ax 

♦ 

♦i 

£ 

n 


Hydraulic  radius 
Eddy  head  loss  term 
Friction  slope 
Slope  of  channel  bed 
Time 

Average  flow  velocity 
Longitudinal  distance 
Vertical  distance 
Flow  depth 

Elevation  of  channel  bed 

Change  in  water-surface  elevation 

Computational  time-step 

Bed  layer  thickness 

Change  in  height  of  natural  levees 

Change  in  height  of  wash  load  deposited  on  floodplain 

Energy  correction  factor 

Momentum  correction  factor 

Density  of  sediment-water  mixture 

Density  of  sediment 

Density  of  water 

Shear  stress 

Derivative  with  respect  to  time 
Derivative  with  respect  to  channel  distance 
Arbitrary  variable 

Variable  evaluated  at  it  net  point  at  the  n**1  time-step 
Width  function 
Integration  variable 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsinile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Johnson,  Billy  H. 

Development  of  a  numerical  modeling  capability  for  the 
computation  of  unsteady  flow  on  the  Ohio  River  and  its 
major  tributaries  /  by  Billy  H.  Johnson  (Hydraulics 
Laboratory,  U.S.  Army  Engineer  Waterways  Experiment 
Station).  --  Vicksburg,  Miss.  :  The  Station  ;  Springfield, 
Va.  ;  available  from  NTIS,  1982. 

154  p.  in  various  pagings  :  ill.  ;  27  cm.  — 

(Technical  report  ;  HL-82-20) 

Cover  title. 

"August  1982." 

Final  report. 

"Prepared  for  U.S.  Army  Engineer  Division,  Ohio  River." 

Bibliography:  p.  46. 

1.  Computer  programs.  2.  FLOWSED  (computer  program). 

3.  Mathematical  models.  4.  Ohio  River.  I.  United 
States.  Army.  Corps  of  Engineers.  Ohio  River  Division. 
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